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ABSTRACT 

Context. Stars are formed in dense molecular clouds. These dense clouds experience radiative feedback from UV photons, and X- 
ray from stars, embedded pre-stellar cores, YSOs, and ultra compact H ll regions. This radiative feedback affects the chemistry and 
thermodynamics of the gas. 

Aims. We aim to probe the chemical and energetic conditions created by radiative feedback through observations of multiple CO, 
HCN and HCO^ transitions. We measure the spatial distribution and excitation of the dense gas (n(H 2 ) > lO'* cm“^) in the core region 
of M17 SW and aim to investigate the influence of UV radiation fields. 

Methods. We used the dual band receiver GREAT on board the SOFIA airborne telescope to obtain a 5'.7x3'.7 map of the J = 
16 —> 15, y = 12 —> 11, and 7 = 11 ^ 10 transitions of '^CO in M17 SW. We compare these maps with corresponding APEX 
and IRAM 30m telescope data for low- and mid-T CO, HCN and HCO^ emission lines, including maps of the HCN 7 = 8 —> 7 and 
HCO^ 7 = 9 —» 8 transitions. The excitation conditions of *^CO, HCO^ and HCN are estimated with a two-phase non-LTE radiative 
transfer model of the line spectral energy distributions (LSEDs) at four selected positions. The energy balance at these positions is 
also studied. 

Results. We obtained extensive LSEDs for the CO, HCN and HCO^ molecules toward M17 SW. These LSEDs can be fit simul¬ 
taneously using the same density and temperature in the two-phase models and to the spectra of all three molecules over a ~ 12 
square arc minute size region of Ml7 SW. Temperatures of up to 240 K are found toward the position of the peak emission of the 
*^CO 7 = 16 —> 15 line. High densities of 10^ cm“^ were found at the position of the peak HCN 7 = 8 —> 7 emission. 

Conclusions. We found HCO^/HCN line ratios larger than unity, which can be explained by a lower excitation temperature of the 
higher-7 HCN lines. The LSED shape, particularly the high-7 tail of the CO lines observed with SOFIA/GREAT, is distinctive for the 
underlying excitation conditions. The cloudlets associated with the cold component of the models are magnetically subcritical and 
supervirial at most of the selected positions. The warm cloudlets instead are all supercritical and also supervirial. The critical magnetic 
field criterion implies that the cold cloudlets at two positions are partially controlled by processes that create and dissipate internal 
motions. Supersonic but sub-Alfvenic velocities in the cold component at most selected positions indicates that internal motions are 
likely MHD waves. Magnetic pressure dominates thermal pressure in both gas components at all selected positions, assuming random 
orientation of the magnetic field. The magnetic pressure of a constant magnetic field throughout all the gas phases can support the total 
internal pressure of the cold components, but it cannot support the internal pressure of the warm components. If the magnetic field 
scales as 6 oc {jjg evolution of the cold cloudlets at two selected positions, and the warm cloudlets at all selected positions, 

will be determined by ambipolar diffusion. 
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1. Introduction 


Inhomogeneous and clumpy clouds, or a partial face-on illu¬ 
mination in star-forming regions like M17 SW, S106, S140, 
the Great Nebula portion of the Orion Molecular Cloud, 
and the NGC 7023 Nebula, produce extended emission of 
atomic lines like [C i] and [C ii], and suppress the stratifi¬ 
cation in [C It], [C i] and CO emission expected from stan¬ 
dard steady-state mode l s of photon-domin a ted regions (PDRsl 


Gerin & Phillins 19981: 

Yamamoto et alJ 12001 : Schneider et al.l 

2002 

20031 iMookeriea 

et alJl2003l: 

Perez-Beaunuits et alJl2010l. 


2012I) . The complex line profiles observed in optically thin 
lines from, e.g., [C i], CS, C'^O, '^CO and their velocity- 
channel maps, are indicative of the clumpy structure of 
molecular clouds and allow for a robust estimation of 


Send ojfprint requests to: J.R Perez-Beaupuits 


their clump mass spectra (e.g. _ 1^^ Lorenl 1^^ ; 


Stutzki & Gustenll99rt Hobsonll99^lRr amer et al 12001 

Perez-Beaupuitsetal. I 2 OI 0 I) . In fact, clumpy PDR models 
can explain extended [C i] ^P\ —> emission and a 

broad LSED for various PD R s and even the whole Milk y 
Wav (IStu tzki & Gtistenl ~990l: Sn aans & van Dishoei^ 1 199 


ICubick et al.ll2008l:l0ssenkopf et alJl2010l). 


Massive star forming regions like the Omega Nebula Ml7, 
with a nearly edge-on view, are ideal sources to study the clumpy 
structure of molecular clouds, as well as the chemical and ther¬ 
modynamic effects of the nearby ioniz ing sources. At a (trigono¬ 
metric paralax) distance of ~ 1.98 kpc (IXu et al.l201 ih . the south¬ 
west region of M17 (Ml7 SW) concentrates molecular material 
in a clumpy structure. The structure of its neutral and molecular 
gas appears to be dominated by magnetic rather than by thermal 
gas p ressure, in contrast to many other PDRs (iPellegrini et al.l 
1200-^ . since magnetic fields as strong as -700 pG have been 
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measured from H i and OH Ze eman observations dBroean et al.l 
[l9^lBrogan & Trolandl2001h . 

Models based on far-IR and sub-m illimeter observations 
dStutzki et al.lll98^ iMeixner et al.lll99^ suggest that the distri¬ 
bution and intensity of the emissions observed in the M17 SW 
complex can be explained with high density (n(H 2 ) ~ 5 x 
10^ cm“^) clumps embedded in a less dense interclump medium 
(n(H 2 ) ~ 3 X 10^ cm“^) surrounded by a diffuse halo (n(H 2 ) ~ 
300 cm"^). 

More recent results from the Stratospheric Observatory for 
Infrared Astronomy (SOFIA) and the German REceiver for As¬ 
tronomy at Terahertz frequencies (GREAT) showed that the line 
profiles and distribution of the [C ii1 158 um emission suppor t 
the clumpy scenario in M17 SW dPerez-Beaupuits et al.ll20l3) 
and that at least 64% of the [C ii] 158 pm emission is not associ¬ 
ated with the star-forming material trace d by the [C i] and C'^O 
lines dPerez-Beaupuits et alj|2013ll2015l) . 

There are more than 100 stars illuminat ing M17 SW, of 
which its central cluster is NGC 6618 (e.g. iLada et alJfToOll : 
iHanson et alJll997 ). The two components o f the massive binary 
CENl dKleinmanr 119731 : IChini et ^Il980h are part of this cen¬ 
tral cluster. This s ource, originally cl assified as a double O or 
early B system by iKleinmarml dl973h . is actually composed of 
two 04 visual binary stars with a separation of ~ 178, named 
CEN la (NE component) and CEN lb (SW component), and 
it appears to be the dominant source of ph oto-ionization in the 
whole M17 region (iHoffmeister et al.ll2008h . 

The components of the CENl binary are also the bright- 
est X-ray sources detected with Chandra in the M17 region 
dBroos et al.ll2007h . These Chandra observations, together with 
a new long (~ 300 ks) exposure of the same field, yielded a 
combined dataset of ~ 2000 X-ray sourc es, with many of them 
having a stellar cou nterpart in IR images. (ITownslev et al.ll2()03l : 
iGetman et al.l 1201^ . Although the brightest X-ray sources are 
CEN la and CEN lb from the central NGC 6618 cluster, other 
stellar concentrations of ~ 40 heavily obscured iE,„edian ^ 2.5 
keV, Ay > 10 mag) X-ray sources are distributed along the 
eastern edge of the M17 SW molecular core. The densest con¬ 
centration of X-ray sources coincides with the well known star¬ 
forming region Ml7-UCl and the concentrati on ends around the 
Kleinmann-Wright Object dBroos et alJl2007l their Eig. 10). The 
M17- UC1 (G15.04-0.68 ), an ultracompact H ii region first stud¬ 
ied by lEelli et aP dl980h . harbors an ionizing source classified 
as a BO-BO.5 star which, together with its southern companion 
(IRS 5S), correspond to a massive Class I source with an ob¬ 
serva ble evolution in the MIR and radio wavebands dChini et alJ 
120001 ) . 

In this study we present maps of multiple high-7 tran¬ 
sitions of '^CO, ^^CO, HCN and HCO^. We analyze the 
underlying excitation conditions from the LSED of some 
of the clumps fitted to their spectral lines. We are partic¬ 
ularly interested in these molecules because their distribu¬ 
tion and excitation is responsive to energetic radiative en¬ 
vironments, which has led to their extensive usage for, 
e.g., disentangling star formation vs. black hole accretion 
and shocks, in the nuclear region of active galaxies (e.g . 


Sternberg et alJll994 KohnoetalJ 19991 Kohno & et al. 2001; 

Kohnol I2003L l2005t 

Usero et al. 120041 

Perez-Beaunuits et al. 


2007 

; Garcfa-Burillo et al.ll2008il: Loenen et al.ll2008; Kri 

ns et al. 

2008 

Perez-Beaunuits et alJ 120091 Rosenberg et al. 

2014; 


Garcfa-Burillo et al.l2()14l:lviti et alj|2014l) 


Given the relat ive youth (< 1 Myr, e.g. iLada et al.l [19911 : 
[Hanson et alJll99'^ of the main ionizing cluster (NGC 6618) of 
M17 SW, and the absence of evolved stars, it is likely that su¬ 


pernovae have not yet occurred. This makes M17 SW an ideal 
place to study the radiative interactions of massive stars with 
their surrounding gas/dust and stellar disks, without the influ¬ 
ence of nearby supernovae. This allows us to study the effects of 
UV, IR, and X-ray radiation fields on the emission of CO, HCN, 
and HCO^ in a well defined environment. 

This study of M17 SW can be considered as one more Galac¬ 
tic template for extra-galactic star forming regions. Therefore, 
we expect that our results will be important for future high 
resolution observations where similar regions in extra-galactic 
sources will be studied in great spatial detail with, e.g., ALMA 
(ISchleicher et al.l[2010l) . 

The organization of this article is as follows. In Sect. 2 we 
describe the observations. The maps of lines observed are pre¬ 
sented in Sect. 3. The modeling and analysis of the ambient con¬ 
ditions are presented in Sect. 4. The conclusions and final re¬ 
marks are presented in Sect. 5. 

2. Observations 

2. 1. The SOFIA/GREAT data 

The new high-7 CO observations were performed with the 
German Receiver for Astrono my at Terahertz Erequencies 
(GREATPl lHevminck et al.l2012h on board the Stratosph eric Ob¬ 
servatory Eor Infrared Astronomy (SOFIA. IVoung et'n]|2012h . 

We used the dual-color spectrometer during the Cycle-1 
flight campaign of 2013 July to simultaneously map a region of 
about 310"x220" (~3.0 pc x 2.1 pc) in the '^CO 7 = 11 10 

at 1381.995105 GHz (216.9 pm) and the ^^CO 7 = 16 ^ 15 
transition at 1841.345506 GHz (162.8 pm) toward M17 SW. 

We also present data for the '^CO 7 = 12 ^ 11 transition 
at 1381.995105 GHz (200.3 pm) that was mapped during the 
Early Science flight on 2011 June. Due to lack of time during 
the observing campaign, the northern map could not be extended 
beyond AS ~ 70". 

The observations were performed in on-the-fly (OTF) total 
power mode. The area mapped consists of six strips, each cov¬ 
ering 224" X 32" (Act x AS with a sampling of 8", half the 
beamwidth at 1.9 THz. Hence, each strip consists of four OTF 
lines containing 28 points each. We integrated Is per dump and 
5s for the off-source reference. 

All our maps are centered on R.A.(J2000)=18*20'"27.6^ 
and Dec(J2000)=-16°12'00'.'9, which corresponds to the star 
SAO 161357. For better system stability and higher observ¬ 
ing efficiency we used a nearby reference position at offset 
(345",-230"). A pointed observation of this reference posi¬ 
tion a gainst the reference (offset: 1040".-535". lMatsuhara et"^ 
Il989h showed that the reference is free of emission. 

Pointing was established with the SOFIA optical guide cam¬ 
eras, and was accurate to a few arc seconds during Cycle- 
1 observations. As backends we use d the Fast Fourier Trans¬ 
form spectrometers dKlein et al.ll2Q12h . which provided 1.5 GHz 
bandwidth with 212 kHz (~0.03 km s ') spectral resolution. 
The calibration of this data in antenna temperature was per- 
formed with the ka librate task from the kosma_software package 
(iGuan et al.ll2012l) . Using the beam efficiencies {rjc) 0.67 for the 
'^CO 7 = 16 ^ 15, 7 = 11 ^ 10 and '^CO 7 = 13 -» 12 lines, 
0.54 for '^CO J - 12 ^ 11, and the forward efficiency (?//) of 


' GREAT is a development by the MPI fur Radioastronomie and the 
KOSMA/ Universitat zu Kdln, in cooperation with the MPI fur Sonnen- 
systemforschung and the DLR Institut fur Planetenforschung. 
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Fig. 1. Colour maps of the integrated intensity of the 7 = 16—» 15,7 = 12—» 11 and 7 = 11 —> 10 transitions of *^CO, and the '^CO 7 = 13 —» 12 
line in M17 SW. The contour levels are the 10%, 25%, 50% (thick line), 75% and 90% of the peak emissions. The central position (Att = 0, 
A6 = 0), marked with a cross, corresponds to the star SAO 161357 at R.A(J2000)=18''20"'27.6‘ and Dec(J2000)=-16°12'00'.'9. The ultracompact 
H II region M17-UC1 and four H 2 O masers dJohnson et all 191^ are marked hy the circle and plus symbols, respectively. The *^CO 7 = 13 —> 12 
map was convolved with a 25" heam to match the resolution of the '^CO 7 = 12 —> 11 and to increase the S/N. All pixels with S/N< 3 (or an 
rms> 0.22 K) were blanked in the *-^CO map. 


0.9'fl we converted all data to main beam brightness tempera¬ 
ture scale, T^h ^rjfX T^/t]c. 

2.2. The APEX data 

We have used the l ower frequency band of the dual channel D SB 
receiver FLASH^ (iHevminck et al.l[2006t iKlein et ^l2014h on 
the Atacama P athfinder Experime nt 12 m submillimeter tele¬ 
scope (APE?^; iGiisten et al.ll2006l) during 2012 November, to 
map the whole M17 region in the 7 = 3 —> 2 and 7 = 4 —> 3 
transitions of '^CO, as well as the '^CO 7 = 3^2 line. The 
observed region covers about 6'.2x7'.2 (4.1 pc x 4.7 pc). In this 
work, however, we present only a smaller fraction matching the 
region mapped with SOFIA/GREAT. 

During 2010 July we performed simultaneous observations 
of HCN and HCO+ 7 = 4 ^ 3 using the 345 GHz band of 
FLASH+ (IKlein et al.l20l4 . as well as the J - 3 ^ 2 transitions 
with the APEX-1 single sideband (SSB) h eterodyne SIS receiver 
(iRisacher et al.ll2006 i lVassilev et al.ll2008l) . 

The regions mapped in HCN and HCO^ cover about 4.5' x 
3.5' (3.0 pc X 2.3 pc) in the 7 = 3 —> 2 transition, and about 
4.5' X 3.2' (3.0 pc X 2.1 pc) in 7 = 4 ^ 3. The maps of HCN 
and HCO^, of both 7 = 4^3 and 7 = 3 —» 2, were done in 
raster mode along R.A. with rows of ~ 270" long (centered on 

^ http://www3.mpifr- bonn.mpg.de/div/submmtech/heterodyne/great/ 
GREAT_calibration.html 

^ APEX is a collaboration between the Max-Planck-Institut fur Ra- 
dioastronomie, the European Southern Observatory, and the Onsala 
Space Observatory 


-95"), with increments of 15" and 10" for the lower and higher 
7-lines, respectively. Subsequent scans in declination with the 
same spacing as in R.A. were done from -120" to h- 90" for the 
7 = 3 —> 2 transition, and from -120" to h- 70" for 7 = 4 ^ 3. 

The total power mode was used for the observations, nod¬ 
ding the antenna prior to each raster to an off-source position 
180" east of the star SAO 161357. The latter is used as reference 
position (Acr -0, A6 = 0) in the maps and throughout the paper, 
with R.A.(J2000)=18''20'”27.6" and Dec(J2000)=-16°12'00'.'9. 
The telescope pointing was checked with observations of con¬ 
tinuum emission from Sgr B2(N) and the pointing accuracy 
was kept below 2" for all the maps. Calibration measurements 
with cold loads were performed every ~10 minutes. The data 
were processed with the APEX real-time calibration software 
(iMuders et al.l 1200^ assuming an image sideband suppression 
of 10 dB for APEX-1. 

For the 7 = 3 —> 2 maps of HCN and HCO^ we used the new 
EFTS with two sub-sections, each 2.5 GHz width, overlapping 1 
GHz in the band center, which gives a total of 4 GHz bandwidth. 
The new FLASH^-345 receiver provides two sidebands sepa¬ 
rately, each 4 GHz wide. Each IF band is then pro cessed with 
two 2 .5 GHz wide backend, with 1 GHz overlap (IKlein et al.l 
I 2 OI 4 I) . We used 32768 channels that gives a spectral resolution 
of about 76 kHz in both 7 = 3 —> 2 and 7 = 4^3 transi¬ 
tions (that is, ~ 8.6 x 10^^ km s ' and ~ 6.6 x 10“^ km s ' for 
the lower and higher 7 line, respectively, at vlsr = 20 kms '). 
The on-source integration time per dump was 10 seconds for the 
HCN and HCO+ 7 = 4^3 maps. For the new FLASH+-345 the 
Tsys was about 180 K, and an average SSB system temperature 
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Fig. 2. Top panels - Maps of the integrated temperature of the HCN and HCO^ 7=1^0 lines in M17 SW. The peak emissions are ~169 Kkm s“' 
and ~152 Kkms^', respectively. The contour levels for all panels in this figure are, from thin to thick, the 10%, 25%, 50%, 75% and 90% of the 
peak emission. Middle panels - Maps of the integrated temperature of HCN and HCO^ / = 3 —> 2 lines, with peak emissions of ~133 Kkm s“' 
and ~'183 Kkms^*, respectively. Bottom panels - Maps of the integrated temperature of the HCN and HCO^ 7 = 4^3 lines in M17 SW. The 
peak emissions are ~'107 Kkms^' and ~187 Kkms^', respectively. The reference position (Act = 0, Ad = 0) marked with a cross, and other 
symbols, are as in Fig.[T] The integrated temperature scale is the same in all maps. 


of 224 K was observed with the APEX-1 receiver for the HCN 
and HCO^ 7 = 3^2 maps. 

We also used the dual color receiver array CHAMP^ on 
APEX to map a region of about 80"x80" (~0.8 pc x 0.8 pc) 
in the HCN 7 = 8 ^ 7 (708.9 GHz) and HCO+ 7 = 9^8 
(802.5 GHz) towards the dense core of M17 SW. The spatial res¬ 
olution of these maps varies between 9.2" for the low frequency 
band, and 8 . 1 " for the high frequency band. 

The coupling efficiency {r]c ~ 0.67) of the new ELASH^- 
345 was estimated from observations at 345 GHz toward Mars, 


during 2010 Juljfl This was extrapolated to 0.659 and 0.656 for 
the HCN 7 = 4^ 3 (354.505 GHz) and HCO+ 7 = 4^3 
(356.734 GHz), respectively. A beam coupling efficiency of 0.72 
was assumed for the A PEX-1 receiver at the frequencies of HCN 
and HCO+ 7 = 3^2 (IVassilev et al.ll200l their Table 2). With 
these beam coupling efficiencies, and a forward efficiency ( 77 /) of 
0.95, we converted all data to main beam brightness temperature 
scale, Tffiij — vjp x 


http://www3.mpifr-bonn.mpg.de/div/submmtech/heterodyne/ 

flashplus/flashmain.html 
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Fig. 3. Line intensity maps of the HCN 7 = 8 —> 7 and HCO'^ 7 = 9^8 transitions in the dense core of M17 SW. The contour levels are the 10%, 
25%, 50% (thick contour), 75% and 90% of the peak emission. The reference position (Aor = 0, AS = 0) marked with a cross, and other symbols, 
are as in Fig. [T]and|2] 



Fig. 4. Left - Velocity-integrated intensity maps of HCN J = 1 —» 0 (grey), HCN 7 = 3 —> 2 (red contour), *^CO J = 11 —» 10 (green contour), 
and *^CO 7 = 16 —> 15 (blue contour). The contour levels (from thin to thick) are the 50%, 75% and 90% of the respective peak emissions. The 
stars on the top-left show the location of O and B ionizing stars (Beetz et al. 1976; Ha nson et al. 1997). Th e reference position (Act = 0, Ad = 0) 
is the same as in Fig[T] The ultracompact H it region M17-UC1 and four H 2 O masers dJohnson et al.lll99^ . are marked hy the (cyan) circle and 
plus symbols, respectively. The horizontal dashed line depict the strip lines shown in Fig. [3 All maps are convolved to a 30" beam, to match the 
resolution of the HCN 7 = 1 —> 0 map. Right - Velocity-integrated intensity maps of HCN 7 = 8^7 (grey), HCO^ 7 = 9 —> 8 (red contour), and 
'^CO 7 = 16 —> 15 (green contour). The blue circles correspond to the heavily obscured (E^cjian > 2.5 keV, Ay > 10 mag) population of X-ray 
sources around the M17-UC1 region (Fig.lO in Broos et al. 2007; coordinates from the VizieR catalogue). In this case, all maps were convolved 
to the beam size (16'.'6) of the *^CO 7 = 16 —> 15 map. 


2.3. The IRAM 30m data 


T he IRAM 30m observa t ions to ward Ml7 SW were described 
in iPerez-Beaupuits et al.l (|2015|) . Here we use the 7 = 1 —> 0 
and 7 = 2 —> 1 maps of '^CO and '^CO reported there. In the 
32 GHz signal bandwidth provided by the I F channels of the 
broadband EMIR receivers dCarter etal .120121) . we also obtained 
data for the 7 = 1 ^ 0 transitions of the HCN, HCO^, as well as 
their isotopologues H'^CN and H'^CO^, in the 3mm band. The 
angular resolutions of the original maps of HCN, HCO^, H'^CN 
and H'^CO^ are 2974, 2972, 3072 and 3071, respectively. 


The reduction of the calibrated data, as well as the maps 
shown throughout the paper, were done using the GILDAI^ 
package CLASS90. All the observed line intensities presented 
throughout this work, as well as the model estimates, are in main 
beam brightness temperature obtained as above, but with a main 
beam coupling efficiencies (rjc) and forward efficiencies {t]/) in¬ 
terpolated for each frequency from the latest reported EMIR ef¬ 
ficiencies We assume Rayleigh-Jeans approximation at all fre¬ 
quencies for consistency. 


5 http://www.iram.fr/IRAMFR/GILDAS 
^ http://www.iram.es/IRAMES/mainWiki/Tram30mEfficiencies 
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Table 1. Line parameters of observed transitions. 


Transition 

Vo" 


E ° 

Wcrtf" 

Telescope/Instnjment 


[GHz] 

["1 

IK] 

[cm“^] 


‘2CO 

7=1^0 

115.271 

22.6 

5.53 

2x10^ 

IRAM 30m/E]VIIR 

7 = 2^1 

230.538 

11.3 

16.60 

7x10’ 

IRAM 30m/EMIR 

J = 3^2 

345.796 

18.9 

33.19 

2x10“ 

APEX/FLASH+ 

7 = 4^3 

461.041 

14.1 

55.32 

4x10“ 

APEX/FLASH+ 

7 = 6^5 

691.473 

9.6 

116.16 

1x10’ 

APEX/CHAMP+ 

J = 7 -i6 

806.652 

8.2 

154.87 

2x10’ 

APEX/CHAMP+ 

7- 11 ^ 10 

1267.014 

24.2 

364.97 

8x10’ 

SOFIA/GREAT 

7 - 12 ^ 11 

1381.995 

22.2 

431.29 

9x10’ 

SOFIA/GREAT 

7 = 13 ^ 12 

1496.923 

20.9 

503.13 

IxlO" 

SOFIA/GREAT 

7= 16^ 15 

1841.346 

16.6 

751.72 

2x10'’ 

SOFIA/GREAT 

'^CO 

7=1^0 

110.201 

23.7 

5.29 

2x10’ 

IRAM 30m/EMIR 

7 = 2^1 

220.399 

11.8 

15.87 

1x10“ 

IRAM 30m/EMIR 

J = 3-^2 

330.588 

19.7 

31.73 

3x10“ 

APEX/FLASH+ 

7 = 6^5 

661.067 

10.0 

111.05 

3x10’ 

APEX/CHAMP+ 

7 = 13 12 

1431.153 

21.4 

481.02 

2x10" 

SOFIA/GREAT 

HCN 

7 = 1^0 

88.632 

29.4 

4.25 

2x10" 

IRAM 30m/EMIR 

7 = 3-i 2 

265.886 

24.9 

25.52 

1x10’ 

APEX/HET230 

7 = 4^3 

354.505 

18.7 

42.53 

3x10’ 

APEX/FLASH 

T 

00 

II 

708.877 

9.2 

153.11 

2x10* 

APEX/CHAMP+ 

h'^cn 

7=1^0 

86.339 

30.2 

4.14 

2x10" 

IRAM 30m/EMIR 

7 = 3^2 

259.012 

25.2 

24.86 

5x10’ 

APEX/HET230 

7 = 4^3 

345.339 

19.2 

41.43 

1x10* 

APEX/FLASH 

HCO+ 

7=1^0 

89.189 

29.2 

4.28 

2x10’ 

IRAM 30m/EMIR 

7 = 3-i 2 

267.558 

24.8 

25.68 

3x10" 

APEX/HET230 

7 = 4^3 

356.734 

18.6 

42.80 

6x10" 

APEX/FLASH 

7 = 9^8 

802.458 

8.1 

192.58 

9x10’ 

APEX/CHAMP+ 

h'^co+ 

7=1^0 

86.754 

30.1 

4.16 

2x10’ 

IRAM 30m/EMIR 

T 

ro 

II 

260.255 

25.1 

24.98 

3x10" 

APEX/HET230 


Notes. Rest frequencies and upper level energies are ado pted from the 
Colog ne Database for Molecular Spectroscopy, CDMS ( iMuller et al.l 
l200.5() and the Lei den Atomic and Molecular Database, LAMDA 
JSchbier et al.ll200^ . ® Beam size (resolution) used to create the re¬ 
spective map. This heam is ~6% larger than the actual HPBW at the 
corresponding frequencies, as defined by the kernel used in the grid- 
ding algorithm of GILDAS/CLASS. Critical densities for temper¬ 
ature ranges 40-300 K (‘^CO, ‘='CO) and 10-30 K (HCN, H'^CN, 
HCO^, H'^CO^). As the main collisional transitions for HCN occur 
with A7 = 2 the traditional two-level formula for the critical density is 
not appropriate here. It generally overestimates the critical density. This 
also applies to the other molecules where AJ = 2 transitions are at least 
comparable. A better estimate for the actual population of the molecule 
is obtained here by adding the coefficients for A7 = 1 and AJ = 2. 


All the spectral lines used and reported in this work are sum¬ 
marized in Table[T] including the upper-level energy and the crit¬ 
ical densities associated with each transition. 


3. Results 

3.1. The high-J line intensity maps 

Figure [T] shows the velocity-integrated (over 10 and 28 kms^') 
temperature (intensity) maps, of the (from top to bottom) J - 
16 ^ 15, 7 = 12 ^ 11, and 7 = 11 —> 10 transitions of '^CO. 
Their peak intensities are about 85 Kkms 316 Kkms ' and 


B 

4<1 




too 0 -too -300 -300 

Aa [ " ] 


Fig. 5. Strip lines of the velocity integrated intensities of all the *^CO 
{top) and the HCN and HCO^ (bottom) lines we have mapped so far, 
at the same declination A6 = 30" (P.A. 90°) across the ionization 
front of M17 SW. The ‘^CO 7 = 7 ^ 6 and 7 = 6 ^ 5 are 
from iPere z-Beaup uits et all ( l2010ll . the '^CO 7 = 13 ^ 12 is from 
iPerez-Beaup uits ej_al. 1 201^jWhile *^CO 7 = 2^1 and 7=1^0 
are from lPdrez^BeaupurtTet all ( 120151) . The rest of the lines are from 
this work. All maps were convolved to a 25" beam to match the reso¬ 
lution of the '^CO 7 = 11 —> 10 map. The X-axis corresponds to the 
actual offset in R.A. of the maps shown in Figs. [T] |2] and The off¬ 
set, Aa = 0" in R.A. correspond to the reference illuminating star SAG 
161357. 


281 Kkms respectively. The '^CO 7 = 12 —> 11, and 7 = 
11^10 follow a similar spatial distribution, and their peaks are 
found at about the same offset position (Aa = -54", A6 = 17"). 
Their emiss ion peaks in b e tween three of the four H 2 O masers 
reported by iJohnson et al.l (1 19981 their Table 9), approximately 
0.3 pc (~ 3073 or 2.5 pixels at P.A. 90°) from the ridge defined 
by the 25% contour line. 

On the other hand, the '^CO 7 = 16 —> 15 emission is much 
fainter and more spatially confined than that observed in the 7 = 
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12^11 and J - 11 —» 10 transitions, and its spatial distribution 
is shifted towards the north-east, closer to the ionization front, 
with respect to the lower-7 lines. In fact, the 7 = 16 ^ 15 
emission peaks very close (~ 12" or 0.1 pc) to the ultra compact 
H II region M17-UC1, which indicates that the '^CO 7 = 16 —> 

15 transition may be excited in warm gas facing the ultracompact 

H II region, not surprising in view of its high upper-level energy 
(£„ 751.72 K). Note that the peak of the '^CO 7=13^12 

closely follows that of the '^CO J - 16 —> 15 emission. 

3.2. The HCN and HCO^ line intensity maps 

Figure|2]shows (from top to bottom) the velocity-integrated (be¬ 
tween 0 and 40 kms ’to cover the hyperhne structure lines of 
HCN) temperature maps of the 7=1—>0, 7 = 3—>2 and 
7 = 4^3 transitions of HCN and HCO^. The peak intensities 
(in Tinb scale) of the 7=1^0 lines of HCN and HCO^ are 
about 169 Kkms ’ and 152 Kkms ’, respectively. The HCN 
andHC0^7 = 3 —» 2 lines have peak intensities ofl33Kkms ’ 
and 183 Kkms ’, and the 7 = 4 —> 3 lines are 107 Kkms^’ 
and 187 Kkms ’, respectively. The HCN 7=1^0 emission 
is brighter than the HCO^ 7=1^0. The opposite is true 
for the higher-7 transitions in agreement with the higher criti¬ 
cal densities of HCN. The HCN and HCO+ 7 = 4^3 lines 
(ricr ~ 2 X 10^ cm“^ and ~ 9 x 10® cm“^ at 100 K, respectively, 
and both with £„ 43 K) are expected to probe much denser and 

colder regions than the (7i,p > 6) CO lines which have £„ 116 K 
(c.f.. Table [T]i. We discuss further the HCO^/HCN line intensity 
ratios in Sec. 13.41 

The HCO^ 7 = 4 —> 3 emission seems more extended than 
the HCN one, particularly towards the northern edge of the cloud 
core, where the ratio is larger. The broader emission of the HCO^ 
7 = 4^3 line w as also observed in the 7 = 3 —» 2 transition 
by iHobsonI dl992l) . who reported a more extended emission of 
HCO^ further into the H ii region. Despite the difference in ex¬ 
tension, the overall spatial distribution (morphology) of the HCN 
and HCO^ emission is similar. 

The HCN 7 = 8 ^ 7 and HCO+ 7 = 9^8 maps 
obtained with CHAMP+ on APEX are shown in Fig. [3 Be¬ 
cause these lines are weak and the beam size of the telescope 
is smaller (with an angular resolution of 9'.'2 and 8'.'1, respec¬ 
tively) we only covered the area around the peak of the HCN 
and HCO^ 7 = 3^2 emission. The peak intensities of the 
HCN 7 = 8 ^ 7 and HCO+ 7 = 9^8 lines are 24 Kkms ’ and 
49 Kkms ’, respectively, and they are located at offset position 
(Aa = -73", A6 - 37"), about 9" (~0.09 pc) from one of the 
H 2 O masers. Because the upper-level energy of HCN 7 = 8 —> 7 
and HCO+ 7 = 9 ^ 8 are £„ 153 K and £„ 193 K, respec¬ 

tively, they may be excited by warmer gas than the lower-7 lines 
of HCN and HCO^, which can explain the second emission peak 
(observed in both lines) found just next to the M17-UC1 region. 

3.3. Comparison of the morphology 

Comparing with the higher-7 ’^CO lines, the left panel of Fig. |4] 
shows the HCN 7 = 3^2, ’^CO 7 = 11 ^ 10 and 7 = 

16 ^ 15 lines overlaid on the HCN 7 = 1 —> 0 map. The 
HCN 7 = 3^2 line follows the distribution of the HCN 7 = 
1^0 line, having their peak intensities at about the same offset 
position. Instead, ’^CO 7 = 11 ^ 10 and ’^CO 7 = 16 ^ 15 
show a clear stratihcation. The peak intensity is shifted towards 
the east relative to the HCN lines (by about 48" or ~0.45 pc 
for ’^CO 7 = 16 ^ 15). The ’^CO 7 = 16 ^ 15 emission is 


very compact and peaks close to the M17-UC1 region. Since the 
7 = 16 ^ 15 and 7 = 11 —> 10 were observed simultaneously 
with the dual band receiver GREAT onboard SOFIA, the shift 
between both lines cannot be due to pointing errors. A difference 
in the excitation conditions is the most likely reason for the shift 
between the peak intensities of these ’^CO lines. 

In the 30" resolution maps, the ’^CO 7 = 16 ^ 15 line peaks 
right in between the M17-UC1 region and two of the four H 2 O 
masers (Fig. 01 left panel). In the higher resolution (1676) map 
shown in the right panel of Fig. 0] the ’^CO 7 = 16 —> 15 emis¬ 
sion actually peaks closer to the UCl region and the eastern most 
H 2 O maser, as well as three of the heavily obscured (E,„edian > 
2.5 keV, Ay > 10 mag) pop ulation of X-ray so urces found 
around the M17-UC1 region bv iBroos et al.l j[2007i their Fig. 10, 
coordinates from the VizieR online catalog. Because of the 
proximity to the ultra compact H 11 region, the H 2 O maser, and 
at least three embedded X-ray sources, the ’^CO 7=16^15 
line is very likely excited by warmer gas than the 7 = 11 —> 10 
transition. This agrees with the fact that the upper-level energy 
{Eup ~ 752 K) of the higher-7 line is about a factor two higher 
than that of the 7 = 11 —» 10 line {E^p ~ 365 K). 

Figure |5] shows the variation of the integrated intensity of 
the ’^CO (top) and the HCN (bottom) lines across the ionization 
front. This corresponds to th e strip line at RA=90° shown by 
(iPerez-Beaupuits et al.ll2010l the ’^CO, ’^CO, 7 = 6^5 and 
7 = 7 —> 6 are taken from that work). The cut also quantihes 
the stratihcation already discussed for Fig. 4. The hrst peak ap¬ 
pears in the high-7 CO lines at -50" while it falls at about -70" 
for all HCN lines. In contrast to ’^CO, HCN shows no shift of 
the peak as a function of the energy level. The peak of the in¬ 
tegrated intensities of the ’^CO lines shows a clear progression 
from deep into the molecular cloud (from about Aa = -120") 
in the 7 = 1 —» 0 and 7 = 2^1 lines towards the molecular 
ridge and the ionization front (Aa between -50" and -60") in 
the 7 = 12 ^ 11 and 7 = 13 ^ 12 lines. The 7 = 16 ^ 15 
emission shows a smooth increment from the molecular cloud 
and reaches a plateau about 10" (~0.096 pc from the peak of the 
7 = 13 ^ 12 line) closer to the ridge at Aa = -40". This incre¬ 
ment in the ’^CO 7 = 16 —> 15 may be indicative of gas getting 
denser and warmer towards the ionization front and close to the 
UC-H II region. On the other hand, the HCN and HCO^ lines 
show a well dehned peak intensity closer to the ionization front, 
between Aa = -60" and -80", coinciding with the peak of the 
’^CO 7 = 7 —» 6 and 7 = 11 —» 10 lines. There is a secondary 
peak deeper into the molecular cloud at Aa between -120" and 
-140", which coincides with the peak intensity of the lower-7 
’^CO lines. 

3.4. The HC0^4HCN line intensity ratios 

The HCO^ 7 = 4 —» 3 map was convolved with the slightly 
larger beam size (17.7") of the HCN 7 = 4^3 line, in order 
to have the same number of pixels per map. The pixel size is 
8.9" X 8.9" which ensures Nyquist sampling of the spectra in 
both (R.A. and Dec.) directions. Likewise, the HCO^ 7 = 3^2 
map was convolved with the beam size (23.6", and spatial reso¬ 
lution of 11.8" X 11.8") of the HCN 7 = 3^2 line before com¬ 
puting the ratio between the respective maps. The HCO^ 7 = 
9^8 was convolved to the larger beam of the HCN 7 = 8^7 
line to match its spatial resolution (9.2"). The same was done for 
the 7=1^0 maps, but convolved to a 30" beam. 


’’ http://vizier.cfa.harvard.edu/viz-bin/VizieR?-source=J/ApJS/ 
169/353 
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Fig. 6. Maps of the HCO^/HCN ratio between the 7 = 1 —> 0 (top-left), 7 = 3 —> 2 (top-right) and 7 = 4 —> 3 (bottom-left) line integrated 
intensities in the region where the lines are brighter than 7% of their peak values. The 7 = 1 —> 0 ratio varies between about 0.4 and 0.9. The 
isolated ratios shown at Ad ~ 70" correspond to significant emission of ^30 Kkms^' in both HCN and HCO^ 7 = 1 —> 0 lines. The 7 = 3 —> 2 
line ratio ranges from ~ 0.8 to ^ 1.7 while the 7 = 4 —> 3 line ratio varies between ^ 1.0 and ~ 2.3. The contour lines correspond to the HCO^ 
7=1—»0, 7 = 3—>2 and 7 = 4 —> 3 maps in the respective line ratios, and are defined as in Fig.|2] The HCO^ 7 = 9 —» 8 to HCN 7 = 8 —» 7 
line ratio is shown in the bottom-right panel, for the region where the lines are brighter than the 15% of their peak values. This ratio varies between 
about unity to ~4.8. The contour lines correspond to the HCN 7 = 8 —> 7 maps and are the same as in Fig.|3 The area of the higher-7 HCN and 
HCO^ ratio map is depicted with a dashed frame in the bottom-left panel. The maps of the higher frequency lines were convolved to the larger 
beam of the lower frequency line before determining the respective ratios. 


In order to reduce uncertainties and avoid misleading values 
in zones with relatively low S/N and no clear detection, we com¬ 
pute the HCO+/HCN 7 = 4^3, 7 = 3^2 and 7=1^0 line 
intensity ratios only in the region where the integrated temper¬ 
ature of both lines (in both species) is brighter than 7% of their 
peak values. That is, between 5cr and 8cr, where cr is the rms of 
the HCN and HCO^ 7 = 4 —» 3 and 7 = 3 —» 2 maps, which 
ranges between 0.2 K and 0.3 K. In the case of the 7 = 1 —> 0 
lines the rms is one order of magnitude lower. Due to the lower 
S/N of the HCN 7 = 8 —> 7 (the fainter line), we used a thresh¬ 
old of 15% for the higher-7 line ratio. This corresponds to about 
6 cr in both HCN 7 = 8 7 and HCO+ 7 = 9 8, with rms 

values of 0.62 K and 1.24 K, respectively. Because the HCN and 
HCO^ maps were observed simultaneously with practically the 
same beam size, our results are not affected by relative pointing 
errors. 

All the HCO^/HCN line intensity ratios are shown in Fig. |6] 
The HCO^/HCN 7=1^0 line ratio is lower than unity in most 
of the region mapped. The fact that the higher-7 HCN/HCO^ 
intensity ratio becomes larger than unity may be due to infrared 


pumping of the HCN 7 = 1 —» 0 line through a vibrational 
transition at 14 pm wavelength. This mechanism can produce an 
enhanced HCN 7 = 1 —> 0 emission in conditions of sub-thermal 
excitation, as discussed in the next section. 

The HCO^/HCN 7 = 3 —> 2 line intensity ratio ranges be¬ 
tween ~ 0.8 and ~ 1.7, while the intensity of the HCO^ 7 = 4^ 
3 line is a factor ~ 1 - 2 brighter than that of the HCN line, as 
shown in Fig.|6] The line profiles of both lines match very well in 
all the positions, which indicates that the difference in intensities 
is not a consequence of different line widths. 

With a higher spatial resolution (972) the line intensity ra¬ 
tio HCO^(9 —> 8)/HCN(8 —> 7) shows the largest values (up 
to ~4.8) of the four ratio maps. The ratio is lowest (between ~1 
and ~2) at the region of the strongest HCN 7 = 8 —> 7 emission 
(as shown by the contours in the bottom-right panel of Fig. 01, 
which coincides with the M17-UC1 region (south-east peak) and 
the peak of the HCN 7 = 8^7 line (north-west peak). The 
relatively bright HCN emission in these two regions can be ex¬ 
plained by dense and warm gas, and by the presence of hot cores 
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(or UC-H II regions) (e g.JPrasad et alJll9^ ilCaselli et al.ll 19931: 


IViti & Williarnsiri99^lRodgers & CharnlevI 2001h. 


Similar high HCO^/HCN 7 = 4^3 line ratios have been 
observed in othe r Galactic star forming regions (e.g., W49A 
iPeng et al.ll20^ . while HCO+/HCN 7 = 4 -» 3 line ratios 
lower than unity have been found from single dish observations 


of active galaxies like NGC 1068 (lPerez;Beau£uits_e^^ 


and in recent ALMA observations (iGarcfa-Burillo et al. 


12009) 


2014 


IViti et all 1201 4 note that for the ALMA observations the au¬ 
thors quote the inv erse HCN/HCO ^ line r atio instead). Accord¬ 
ing to models by iMeiierink et alJ (l2007l their Fig. 14), how- 
HCO^/HCN 7 = 4^3 line ratios lower than unity 


ever. 


are signatures of a high density (> 10“* cm“^) PDR environ¬ 
ment rat her than an X-ray dominat ed region (XDR), as sug¬ 
gested bv iGarcfa-Burillo et al.l (l2014l) . But interpretations of the 
HCO^/HCN line ratios in extragalactic environments may not 
be as straight forward as believed in the past, since these species 
are now known to be s trongly time-depe ndent in both dense ga s 
(e.g. lBavet et al.l2008h and i n XDR s (e.g. lMeherink et al.l2013h . 
as pointed out bv IViti et al.l (l2014t) . Hence, the HCO^/HCN line 
ratios observed in galaxies may reflect departures from chemical 
equilibrium rather than differences in the underlying excitation 
environment. 


3.4.1. Infrared pumping of the dense molecular gas? 

In environments with high enough UV/X-ray surface brightness 
to heat the dust to several hundred K, like in AGNs, bright 
starbursts or massive star-forming regions in the Galaxy and 
the circumnuclear disk (CND) around the massive black hole 
at the center of our Galaxy, strong HCN, HNC, and HCO^ 
emission in the sub-millimeter and millimeter range may also 
be explained by t he infrared radiative pumping scenario (e.g. 


Aalto et al. 19951 Christo 

pheretal. 2005; Garcfa-Burillo et alJ 

20061 iGuelin et al. 12007; 

Aalto et al.| 2007b| al). since hot dust 


produces strong mid-infrared continuum emission in the 10-30 
pm wavelength range. 


If the infrared pumping scenario is at work, absorption fea¬ 
tures must be detected in infrared spectra at 12.1 pm, 14.0 pm, 
and 21.7 pm, which are the wavelengths of the ground vibra¬ 
tional states of HCO^, HCN, and HNC, respectively, connected 
by the transitions V 2 = 1 of the lowest excited bending states. At 
the moment we do not have high resolution IR spectra to check if 
this is the case in all the mapped region of M17 SW. On the other 
hand, the HCN 7 = 4 —» 3, V 2 = 1 transition (356.256 GHz) was 
detected with a peak intensity of ~50 mK towards the CND of 
the Milky Way (iMills et al.ll2013l) . indicating radiative pumping 
of HCN at 14.0 pm. However, we did not detect this vibrationally 
excited HCN line, which lies close to our HCO^ spectra, at an 
rms level of ~ 1 mK. Besides, the fact that the HCO^/HCN line 
ratio grows monotonically with 7-level can be considered an ar¬ 
gument against 14 pm pumping. Therefore, we think the HCN 
and HCO^ lines in M17 SW are unlikely to be affected by in¬ 
frared pumping. 


3.4.2. Other causes of high HCOVHCN line ratios 

Due to the smaller physical scales and the different metallicity, 
the time-dependency of the abundance of HCN compared to that 
of HCO^ should not be an issue in Galactic molecular clouds. 
However, the bulk HCO^ emission may still arise from gas that 
does not co-exist (in terms of ambient conditions or, equiva¬ 
lently, from a different layer within a clump) with the gas host¬ 


ing the HCN, as previously suggested bv IPerez-Beaunuits et al.l 
(l2()09l their Apendix C.3). 

Stronger HCO^ emission could also be the consequence 
of the higher ionization degree in X-ray dominated regions 
(XDRs), which leads to an enhanced HCO^ formati on rate (e.g. 
iLepp & Dalgarnolll996l: IMeiierink & Spaansll2005l) . Using the 
XDR models bv IMeiierink & Spaan^2005 ) we estimate that in 
order to drive an XDR, an X-ray source (or a cluster of sources) 
with a (combined) impinging luminosity of at least 10^^ erg s“' 
would be required to be within a few arcsecs (< 0.03 pc in 
M17 SW) from the region of the HCO^ 7 = 4^3 peak emis¬ 
sion, considering that the X-ray flux decreases with the square 
of the distance from the source. 

Based on the luminosities estimated from thermal plasma 
(546 sources) and power law fits (52 sou rces) of the photorn etri- 
cally selected ChandrafAClS sources bv iBroos et al.l (l2007h . we 
estimate that the combined X-ray luminosity (integrated between 
0.5 keV and 50 keV, and projected on the sky at the position 
of the peak HCO^ emission) is about three orders of magnitude 
lower than needed to drive an XDR. This result rules out an XDR 
scenario to explain the relatively strong HCO^ emission. How¬ 
ever, our estimates are based on an X-ray SED fit, extrapolated 
up to 50 keV with no actual observations above 10 keV (the up¬ 
per band of ChondralACIS). Higher energy (>10 keV) photons 
could make a larger contribution to the X-ray flux, particularly 
from sources that have significant power-law tails. However, un¬ 
less such X-ray sources would be located sufficiently close to 
(or within) the region with bright HCO^ emission, they proba¬ 
bly could not increase the required X-ray flux by three orders 
of magnitude. In fact, no X-ray source was detected by Chan¬ 
dra/ACIS within a radius of ~10" (~0.09 pc) around the peak 
HCO^ emission. This either rules out the existence of any X-ray 
source in that region or implies that all X-ray photons with en¬ 
ergy <10 keV are heavily a bsorbed by the large col umn density 


_ large 

of the gas (~8xl0^^ cm~^. IStutzki & Gustenlll9^ . Therefore, 
future observations sensitive to higher energy {kT > 10 keV) 
photons that could escape the large column of gas, are required to 
unambiguously discard any heavily obscured X-ray source that 
may exist in the dense core of M17 SW. 

Besides X-rays there are other possible mechanisms that can 
also lead to a relatively bright HCO^ emission. Indeed, if the in¬ 
tensity of emission from HCO^ is as strong as that from HCN (or 
stronger) it may be due to relatively high kinetic temperatures, 
strong UV radiation fields, and relatively lower de nsities of the 
gas from where the HCO^ emission emerges (e.g. Fuente_e^^ 


1993; Onn^t^ , 12221 ’ rBrouillet et al.l 200 ^ IChristopher et ^ 

TZb- 


2005:Eiang et al. 120071: IMeiierink et al.ll200 

Since HCN and HCO^ are expected to co-exist at similar 
depths in a cloud, the kinetic temperature of their surrounding 
gas must also be similar. But the critical density of the HCO^ 
lines is one order of magnitude lower than that of HCN (cf., Ta- 
ble[TJ, which makes HCO^ more easily excited (collisionally) in 
single-phase molecular gas. Because the HCO^/HCN line ratio 
reflects the effects of the combination of abundance and excita¬ 
tion temperature, the increasing ratio with 7-line can be indica¬ 
tive of a lower excitation temperature in the higher-7 transitions 
of HCN. We explore this alternative through excitation models 
in SectlH 


4. Analysis of the LSEDs 

We first computed the average spectra of a 200 arcsec^ region to 
mimic the results obtained when observing extragalactic sources 
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Fig. 7. Velocity-integrated intensity map of '^CO J - 3 2, convolved 

a 25" beam to match the resolution of the '^CO 7 = 11 —> 10 map. The 
contour levels are the 10%, 25%, 50% (thick line), 75% and 90% of the 
peak emission (770 Kkms^*). The central position (Aa = 0, Ad = 0) 
marked with a white cross, is the same as in Fig. [T| The ultracom pact 
H II region M17-UC1 and four H 2 O masers i.Iohnson et al.lll99^ are 
marked by the filled circle and plus symbols, respectively. The positions 
of the peak intensities of HCN 7 = 8 —> 7 and *^CO 7 = 16 —> 15, as 
well as the offset positions at (-60",-30") and (-130",-1-30") analyzed 
in Sec. 14.41 are indicated with a triangle, a diamond, a pentagon and a 
star, respectively. 
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Fig. 8. Average spectra of *^CO (top), HCN (middle - indicating the 
hyperfine structure lines of HCN 7 = 1 —> 0 ) and HCO^ (bottom), 
obtained from the 200"-beam area (dashed circle in Fig. |7]( toward 
M17SW, centered at offset position (-100'',0"). This corres pond to a 
regio n of ~1.92 pc in diameter at the distance of 1.98 kpc dXu et all 

HSl!). 


with ALMA. These average spectra can be fit using a single 
Gaussian component, as shown in Fig. [8] 

Then we analyze the spectra from the 25" resolution maps at 
four selected positions in M17 SW. The profiles of the CO lines 
show a very rich structure due to the dumpiness of the source 
(several velocity components along the line of sight) and to opti¬ 
cal depth effects (self-absorption) of the lower-7 transitions (cf.. 
Fig. [Toll. While the lower- and mid-7 '^CO lines can be fit with 
up to five Gaussian components, the higher-7 (Jup >11) can be 
fit with just one or two Gaussian components. In order to en¬ 
sure we are comparing the line fluxes corresponding to the same 
velocity component, we chose the Gaussian component of the 
lower- and mid-7 lines that is the closest related to the central 
velocity and line width obtained for the single Gaussian used to 
fit the '2cO 7=16^15 line. 

The profiles of the HCN and HCO^ lines also show 
structures that could be associated with different overlapping 
cloudlets along the line of sight. However, contrarily to the '^CO 
lines, the H'^CN and H'^CO"^ lines do not show the same struc¬ 
ture as the main isotope lines, (cf.. Figs. [TT] and [T2t . With the 
exception of the 7=1^0 transitions, the higher-7 '^C bear¬ 
ing lines seem to be made of a single component that can be 
fit with one Gaussian. This is indicative of self-absorption (or 
optical depth effects) in the '^C lines. Therefore, instead of us¬ 
ing several Gaussian components in the HCN and HCO^ lines, 
we fit a single component trying to reproduce the missing flux 
by using just the apparently less self-absorbed emission in the 
line wings and masking the line centers for the Gaussian fit. The 
fluxes obtained in this way are significantly larger (between 13% 
and 95%, depending on the line and the offset position) for the 
HCN lines, compared to the flux obtained by using two or more 
Gaussian components. The HCO^ lines show a lower degree of 
self-absorption. 

After obtaining the velocity-integrated intensities from the 
Gaussian fit, we convert the intensities to fluxes in units of 

erg s“' cm“^ sr“' as J T,„i,dv, where k is the 

Boltzmann’s constant, vq is the rest frequency of the transition, 
c is the speed of light and the factor 10^ is to convert from 
kms * to cm s '. This flux is then converted to units of W m^ 
by Fyg - 10“’F'^Q25" 10"' with G 25 " the 25" beam solid angle 


(resolution) used in our maps, the factor 10"' is to convert from 
cm^ to m^, and the factor 10“^ to convert from erg s“' to Watts. 

We model the observed flux es using the modified RADEX0 
code (Ivan der Tak et al.l 1200']^ that uses a background ra- 


;s a ba c kgi 

diation field, as done by iPoelman & Snaansl (120051) and 
jPerez-Beaupuits et al.l (1200 9 ), and conside ring a dust tempera- 
ture of 50 K found bv iMeixner et al.l (Il992h toward M17 SW. 

We chose this modified background radiation field because 
interaction with far- and mid-infrared radiation (mainly from 
dust emission in circumstellar material or in star- forming re¬ 
gions) can be important for molecules with widely spaced rota¬ 
tional energy levels (e.g., the lighter hydrides OH, H 2 O, H 30 ^ 
and NH 2 ), and particularly for the higher-7 levels of heavy 
molecular rotors such as CO, CS, HCN, HCO^ and H 2 CO, that 
are observable with Herschel and SOFIA in the far- and mid- 
infrared regime. However, the actual effect of the dust IR emis¬ 
sion on the redistribution among rotational levels of a molecule 
depends on the local ambient conditions of the emitting gas. That 
is because at high densities (or temperatures) collisions are ex¬ 
pected to dominate the excitation of the mid- and high-7 levels 
of molecules, such as CO, while at lower densities (or temper¬ 
atures) radiative excitation, as well as spontaneous decay from 


http://www.sron.rug.nl/~vdtak/radex/index.shtml 
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higher-7 levels, are expected to be the dominant component driv¬ 
ing the redistribution of the level populations. Since we do not 
know a priori what the ambient conditions are at the selected 
positions in Ml7 SW, we include the modified background radi¬ 
ation field for completeness. 

The RADEX code computes the local escape probability as¬ 
suming a uniform temperature and density of the collision part¬ 
ners. For the case of M17 SW the escape probability can be com¬ 
puted using two independent methods implemented in RADEX; 
the large velocity gradient (LVG) formalism and assuming an 
homogeneous sphere geometry. We tested both methods against 
each other with a test fit of the LSEDs using low enough densi¬ 
ties and column densities to ensure convergence. The results ob¬ 
tained with both methods were not significantly different for the 
CO LSEDs; the rms (or;^f^) values obtained from the observed 
fluxes and those predicted by the two formalisms differ by just 
~2% between the methods, with the largest differences pertain¬ 
ing for Jup > 10 CO lines. For the HCN lines, however, the LVG 
method gives about 34% higher intensities, while for the HCO^ 
lines the LVG method gives between 15% and 20% higher inten¬ 
sities. In order to obtain similar HCN and HCO^ line intensities 
as with the LVG method the uniform sphere formalism requires 
larger column densities and slightly larger filling factors. When 
using relatively large densities (>4x10^ cm“^) and column den¬ 
sities (>2xl0'^ cm“^) for HCN the LVG method did not con¬ 
verge. Since these parameters ranges turned out to be relevant for 
a good fit of the LSEDs, we used the uniform sphere geometry 
method in our models. The physical conditions were modeled 
using the collisiona l data available in the LAMDA0 database 
(ISchoier et al.ll20p5h. The co llisional rate coefficients for '^CO 
are adopted from Yang et akl (120101) . 


4.1. The Spectral Line Energy Distribution 

We had to fit the full LSEDs of the '^CO, HCN, and HCO+ lines 
by assuming a cold and a warm component. A single component 
does not fit all the fluxes and three components require more 
free parameters than the available number of observations (par¬ 
ticularly for the HCN and HCO^ lines). The model we use is 
described by: 


Etot(y) = <l> coldFcold "I" ^warmFwarm (1) 

where O; are the beam area filling factors and L, are the es¬ 
timated fluxes for each component in units of W m“^. The 
estimated fluxes of the main isotopologues (i.e., ^^C-bearing 
molecules) are a function of four parameters per component: 
the beam area filling factor <1), the density of the collision part¬ 
ner n(H 2 ) (cm“^), the kinetic temperature of the gas (K), and 
the column density per line width NfAV (cm^^km^* s) of the 
molecule in study. We used the average FWHM obtained from 
the Gaussian fit of the Jup > 2 lines for the line widths AV of 
a given species (the Jup < 2 transitions were not considered be¬ 
cause their emission was collected with a beam size larger than 

25"). 

The two component fit of the LSED req uires eight free 
parameters. We use the simplex method (e.g., iNelder & MeadI 
Il965t iKolda et^l2003h to minimize the error between the ob¬ 
served and estimated fluxes, using sensible initial values and 
constraints of the input parameters as described below. The un¬ 
certainty in the parameters is obtained by allowing solutions of 
the LSED fit in a 3cr range around the estimated fluxes, where cr 

® http://www.strw.leidenuniv.nl/~moldata/ 


Table 2. Gaussian fit parameters for the average spectra. 


Transition 

f T,„idV 

rppeak 
^ mb 

Vo 

AV (FWHM) 



[Kkms-'l 

[K] 

[km s“' ] 

[kms^'l 

O 

U 

y = 1 ^ 

0 

264.0± 5.2 

30.66±0.95 

19.89± 0.08 

8.09± 0.19 

J = 2-> 

1 

282.9± 1.0 

31.95±0.18 

19.88± 0.02 

8.32± 0.04 

7 = 3^ 

2 

360.4± 0.0 

39.20±0.01 

20.21±0.00 

8.64± 0.00 

7 = 4-> 

3 

328.0± 0.1 

36.62±0.01 

19.86± 0.00 

8.41±0.00 

J = 6~> 

5 

342.7± 0.0 

42.43±0.00 

20.73± 0.00 

7.59± 0.00 

7 = 7-^ 

6 

354.3± 1.0 

45.84±0.20 

20.75± 0.01 

7.26± 0.02 

7=11- 

^ 10 

72.0± 0.5 

12.30±0.15 

20.63± 0.02 

5.50± 0.05 

7= 12- 

A 11 

44.9± 0.7 

8.28±0.22 

19.87± 0.04 

5.10± 0.10 

7= 13 - 

A 12 

34.7± 1.2 

6.69±0.38 

20.65± 0.08 

4.87± 0.21 

7= 16- 

^ 15 

5.9± 0.4 

1.31±0.16 

20.10± 0.16 

4.28± 0.41 

'^co 

7=1^ 

0 

63.0± 0.7 

10.88±0.19 

19.74± 0.03 

5.44± 0.07 

7 = 2^ 

1 

169.0± 0.0 

26.91±0.01 

19.96± 0.00 

5.90± 0.00 

J -3 -* 

2 

157.6± 0.4 

24.35±0.09 

20.09± 0.01 

6.08± 0.02 

7 = 6 —> 

5 

124.6± 0.0 

21.67±0.00 

20.57± 0.00 

5.40± 0.00 

7= 13 - 

^ 12 

1.9± 0.5 

0.28±0.13 

20.39± 0.91 

6.62± 2.55 

HCN 

7=1^ 

0 

40.1± 1.8 

6.51±0.29 

19.15±0.20 

5.80± 0.00 

J=3-> 

2 

28.3± 0.0 

3.65±0.00 

19.84± 0.00 

7.29± 0.00 

7 = 4^ 

3 

20.9± 0.0 

3.00±0.01 

19.98± 0.01 

6.57± 0.02 

7 = 8 —> 

7 

1.4± 0.1 

0.29±0.04 

19.84± 0.25 

4.77± 0.53 

h'^cn 

7=1^ 

0 

5.7± 0.1 

0.60±0.02 

19.61±0.09 

9.01±0.23 

7 = 3-^ 

2 

2.1±0.0 

0.39±0.03 

19.44± 0.11 

5.06± 0.26 

7 = 4-^ 

3 

1.0± 0.0 

0.23±0.01 

19.56± 0.06 

4.47± 0.14 

HCO+ 

7=1^ 

0 

40.7± 0.1 

6.28±0.02 

19.43± 0.00 

6.09± 0.01 

7= 3 -» 

2 

33.2± 0.0 

5.29±0.02 

20.76± 0.00 

5.91±0.02 

7 = 4-^ 

3 

34.9± 0.0 

5.54±0.01 

20.69± 0.00 

5.92± 0.01 

7 = 9-> 

8 

3.3± 0.4 

0.71±0.16 

19.88± 0.25 

4.40± 0.79 

H'3C0+ 

7=1^ 

0 

3.7± 0.0 

0.80±0.02 

19.51±0.03 

4.37± 0.07 

7 = 3-> 

2 

3.4± 0.1 

0.75±0.03 

20.19± 0.06 

4.32± 0.15 


includes the standard deviation of the fluxes obtained from the 
Gaussian fit and 20% of calibration uncertainties. 

We also included all the available '^CO fluxes to constrain 
mainly the column density. We use the same excitation condi¬ 
tions (density and temperature) as for '^CO. The same was done 
for the H'^CN and H'^CO^ lines. In order to reduce the number 
of free parameters, we use the same filling factors for the as 
for the bearing lines in the cold and warm component. This 
is a reasonable assumption since the lines, although usually 
fainter and less extended when shown in linear scale, appear as 
extended as the lines when shown in logarithmic scale. 

Because for HCN and HCO^ we have fewer transitions than 
for CO, we first fit the CO LSED and then used the same density 
and temperature found for the cold and warm components from 
CO for fitting the HCN and HCO+ LSEDs. 

4.2. Constraints on the model parameters 

Two of the most critical parameters of the models are the column 
densities and the beam area filling factors. The filling factors 
were discussed above, while the column densities depend on the 
observed fluxes and the excitation conditions (temperature and 
density) of the models. 

Another important parameter of the models is the carbon- 
12 to carbon-13 isotope ratio, which couples the column den¬ 
sities of the '^C and '^C lines since we assume A(*^CO) = 
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Fig. 9. Two component fit of the line spectral energy distribution of the CO, HCN, and HCO^ species, for the average spectra toward M17 SW. The 
fit of the bearing isotopologues are shown in the insets. The cold and warm components are shown in dotted and dashed lines, respectively. 
The error bars correspond to the 3cr range used to estimate the standard deviation of the excitation parameters, where cr is the flux uncertainties 
obtained from the Gaussian fit. 


Table 3. LSED fit parameters" for the average spectra. 


4.3. Excitation of the Average Emission in M17 Sl/l/ 


Parameter CO HCN HCO+ 


®»m(‘"C) 
ncoldOil) [cm“^] 


Tcold IK] 

Ncold Icm“^] 

«»arm(H2) [cm“^] 
Twarm [K] 

^warm [cm“-] 


1.00 ±0.13 
4.80 ± 0.25 
42.00 ±3.71 
19.50 ±0.48 

0.10 ± 0.01 
6.00 ± 0.60 
135.00 ± 8.46 
18.10 ±0.36 

1.00 ± 0.12 
0.10 ± 0.00 


0.60 ± 0.07 
4.80 ± 0.49 
42.00 ± 3.91 
15.40± 1.19 

0.10 ± 0.01 
6.00 ± 0.77 
135.00 ± 12.09 
14.70 ± 0.94 

0.60 ± 0.06 
0.10 ± 0.01 


0.60 ± 0.07 
4.80 ±0.54 
42.00 ±4.38 
14.20 ±0.88 

0.10 ± 0.01 
6.00 ± 0.48 
135.00 ± 15.51 
14.40 ± 0.50 

0.60 ± 0.05 
0.10 ± 0.01 


[kms-1] 8.00 6.00 6.00 

[kms-1] 6.00 6.00 4.00 


Notes. The density and column density values are given in /ogio scale. 


Ai(i2C0)/'Ri2/i3, with 'Ri 2/13=[‘^C]/['3C]. The isotope ratio 
7 ?i2/i 3 is known to vary from source to source in the Milky 
Way, depending on the distance to the Galactic center (e.g. 


Henkel et alJll982l 1 19851 : iLanger & Penzias 


2005 !) . At a distance of 1.98 kpc JXu et all 


19901: iMilam et alJ 


201 Ih . M17 SW is 


expected to have a CO isotope ratio of *R i 2 /i 3 = 51.6 + 7.9, 
according to Eq. (3) bv iMilam et alJ (l2005h . In our models we 
use I??i 2 /i 3 = 50. For HCN and HCO^ the same isotope ra¬ 
tio as for CO can be considered. However, there is theoretical 
and observational evidence that the actual HCN isotope ratio 
is closer to that obtained from observations of H 2 CO, which 
provides an upper limit to 7?i2/i3 while CO provides a lower 
limit, and is otope ratios derived from HCO^ lead to intermediate 
values (e.g. lLanger et al.lll984t fHenkel et alJll985l: IMilam et ^ 
I 2 OO 5 I 1 . Thus, from Eqs. (5) and f41 bv IMilam et alj l 2005 1. the 
isotope ratios for HCO^ and HCN in M17 SW should be ~56 
and ~63, respectively. Using these higher 7?i2/i3 values, how¬ 
ever, we could not fit the HCN and HCO^ LSEDs simultane¬ 
ously with the H^^CN and H'^CO^ lines. We were able to fit 
the HCN and HCO^ LSEDs using the same isotope ratio of 50 
as for the CO lines, although smaller 7?i2/i3 of 30 and even 20, 
would lead to a better fit of the rare HCN and HCO^ isotopo¬ 
logues. T his would be in agreement w ith the more recent results 
found bv iRollig & Ossenkopfl d2013h for HCO^, but not for the 
HCN isotope ratio. The discrepancy we found for HCN is likely 
due to optically thick lines which are also heavily affected by 
self-absorption, as mentioned in Sect. 0] 


We computed the average spectrum in a region associated with a 
200 " beam size, centered at offset position (-100",0"), in order 
to estimate the excitation derived from ALMA observations of 
nearby Galaxies. In M17 SW the 200" beam covers a spatial 
scale of ~2 pc. This correspond to the spatial scale that would be 
resolved by ALMA with the finest achievable angular resolution 
(0.03" with the capabilities available in Cycle 3) of the bands 6 
(230 GHz) and 10 (870 GHz) towards a galaxy like NGC 1068 
at a distance of ~ 14.4 Mpc. Coarser angular resolution would be 
sufficient to resolve the same spatial scale in closer galaxies like 
NGC 253 (at ~3.5 Mpc). 

Figure I 2 ] shows the ^^CO J - 3 ^ 2 map convolved with a 
25" beam, and the dashed circle depicts the area where the av¬ 
erage emission in all the maps was estimated from. The average 
spectra of the ^^CO, HCN and HCO^ lines are shown in Fig. 0 
We computed the intensities integrating the spectra over the 5- 
35 km s * velocity range. We also fit one Gaussian component to 
estimate the (FWHM) line width. The velocity-integrated inten¬ 
sities and the total intensity obtained from the Gaussian fit are 
similar within a few percent. 

The best two components model fit of Eq. [T] is shown in 
Fig. |9] for I^CO, HCN and HCO+. The LSED fit of the isotopo¬ 
logues, '^CO, H^^CN and H'^CO^, is shown in the insets. The 
error bars correspond to the 3cr range around the fluxes used to 
estimate the standard deviation of the parameters, where cr is 
the uncertainty of the estimated fluxes obtained from the Gaus¬ 
sian fit. Even though the fluxes of the HCN fine structure lines 
are not considered in the Gaussian fit, the flux estimated for the 
HCN 7 = 1 —> 0 line is similar than predicted by the model 
(Fig .|3 middle panel). 

The parameters of the LSED models are summarized in Ta¬ 
ble [H Since the CO emission is more extended than the HCN 
and HCO^, the beam area filling factors of both CO components 
need to be larger as well. The same excitation conditions (density 
and temperature) as found for CO, were sufficient to reproduce 
the HCN and HCO^ fluxes, considering the same filling factors 
for these two species but with 37% and 84% lower column den¬ 
sities for the cold and warm component of HCO^, respectively. 
At first look, this might be counter-intuitive with the fact that the 
emission from the higher-7 (Jup > 3) transitions of HCO^ are 
brighter and seem to be more extended than those of HCN (cf.. 
Figs. |2] and |3]l. But this can be explained because HCO^ is easier 
to excite than HCN, given their different critical densities (see 
Table[T]l, and the larger optical depths of the HCN lines. 
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Fig. 10. Spectra of all '^CO {left panel) and '^CO {right panel) lines observed at the four selected offset (within ±3") positions shown in Fig.|7] 
All the spectra, except the 7 = 1 —> 0 and 7 = 2 —> 1 have been convolved with a 25" beam, corresponding to the '^CO 7 = 11 —> 12 map. 


The densities and temperatures found for the cold compo¬ 
nent of CO are compa rable with those found in starburst galax¬ 
ies like, e.g. NGC 253 dRosenberg et alJl20l4 their Tables 2^), 
but the warm component requires excitation conditions with at 
least one order of magnitude higher density and column density. 
Similar discrepancies are found when comparing with the exci- 
tation conditions estim ated for Seyfert galaxies like NGC 1068 
(ISninoglio et al.l 1201 2l their Table 3). These differences may 
be due to an underestimated beam dilution effect in the ex- 
tragalactic observations, as well as the fact that in their larger 
beams they collect emission emerging from a variety of molec¬ 
ular clouds with different sources of heating. Another factor af¬ 
fecting the determination of the excitation conditions in extra- 
galactic observations is the dust extinction that seems to affect 
the high-7 CO lines at the large (>10^"^ cm~^) column den- 
sities observed towards their centers (e.g. iPineda et al.l 120101 : 


lEtxaluze et al.l 120131) . which was not taken into account in the 
extragalactic studies. 

We do not correct the higher-7 CO lines for dust extinction 
either since the dust extinction in Ml7 SW is expected to be less 
significant than towards circumnuclear regions, given that the to¬ 
tal column densities estimated towards M17 SW are between one 
and two orde rs of magnitude smaller than towards our Gal actic 
Center (e.g., IStutzki & Giistenl [19901: iMeixner et al.lll99^ and 
the centers of other galaxies. 

4.4. Excitation at key positions in M17 Sl/V 

Using the 25" resolution maps we estimated the fluxes from the 
Gaussian fit of the line wings, as described in the third paragraph 
of Sect.|4] Then we fit the USED of the three molecular species at 
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Fig. 11. Spectra of all HCN (leftpanel) and H'^CN (right panel) lines observed at the four selected offset (within ±3") positions shown in Fig. [7] 
All the spectra, except the 7 = 1 —> 0 and 7 = 2 —> 1 have been convolved with a 25" beam, corresponding to the *^CO 7 = 11 —> 12 map. 


four particular locations toward Ml7 SW to study the variation 
in excitation conditions. 

We chose the peak of the HCN 7 = 8^9 emission at 
about the offset position (-65",+31") which is expected to be 
dominated by dense gas, the peak of the 7=16^15 

emission at about (-40",+18"), close to the UC H ii region, the 
emission at about (-60",-30"), deeper into the PDR with an 
expected fair mixture of excitation conditions, and the emission 
at about (-130",+30") west from the northern concentration of 
H 2 O masers, corresponding to the shoulder or secondary peak of 
the HCN and HCO^ strip lines of Fig. |5] We refer to these posi¬ 
tions as: HCN-peak, CO-peak, PDR, and M17-W, respectively, 
and indicated in Fig. |7] 

The spectra of all the observed lines at these positions are 
shown in Figs. [TOlfT^ The fluxes obtained from the CO, HCN 


and HCO^ spectra at the M17-W position have the largest uncer¬ 
tainties because the line profiles show the highest level of sub¬ 
structures and self-absorption features. The H^^CN lines at this 
position have very low S/N, so we consider them as upper limits. 
The best fit models for the LSED at these positions are shown in 
Figs. [ 13 ] and the model parameters are summarized in Tables]?] 

SISIandiTj 

The '^CO LSED at the HCN-peak shows an increasing trend 
in flux with 7, with an apparent turn over at the 7 = 14 ^ 13 
transition (which we do not have). The '^CO LSED of the CO- 
peak emission seems flat at the high-7 transitions, and the '^CO 
LSED at the PDR position, shows a decreasing trend in flux from 
the 7 = 11 —> 10 transition. On the other hand, the LSED of the 
M17-W position show a flat distribution of fluxes between the 
J - 7 ^ 6 and J - 13 ^ 12 transitions, and a sharp (one order 
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Fig. 12. Spectra of all HCO^ (leftpanel) and H'^CO^ (rightpanel) lines observed at the four selected offset (within ±3") positions shown in Fig.|7] 
All the spectra, except the 7 = 1 —> 0 and 7 = 2 —> 1 have been convolved with a 25" beam, corresponding to the *^CO 7 = 11 —> 12 map. 


of magnitude) decrease at the 7=16^15. Observations of the 
7 = 9 —> 8 transition with APEX would be needed to confirm 
the flat distribution of fluxes at the mid- and high-7 transitions 
of at the M17-W position. These different LSED shapes, 
specially at the high-7 CO lines obtained with SOEIA/GREAT, 
are indicative of the distinctive excitation conditions dominating 
mainly the warm component of the two-phase model. 

The highest gas temperature of all selected positions was 
that of the warm component of the CO-peak, ~240 K, which is 
in agreement with temperatures previously estimate d toward the 
PDR i nterface from high- and mid-7 ^^CO lines (e.g. lHarris et al.l 
Il987t IStutzki et al.lll98^ . Although the previous studies used 
larger beam sizes, probably collecting emission from colder gas, 
their estimates lead to a range of likely temperatures between 
100 K and 500 K, and densities of few times 10"^ cm“^. In¬ 
stead the warm component of our two-phase models require be¬ 


tween one and two orders of magnitude higher densities. This 
discrepancy in densities is explained because the previous esti¬ 
mates were based on line ratios of the '^CO 7 = 7 —> 6 and 
7 = 14 —> 13 transitions, and line ratios of only two transitions 
are subject to the know dichotomy between density and temper¬ 
ature of the gas. Our estimates, instead, are constrained in great 
degree by the mid- and high-7 '^CO, H^^CN, and H'^CO^ iso- 
topologue lines, which require high (> 10^ cm“^) densities to 
reproduce the observations. Lower densities (< 10^ cm"^) and 
higher temperatures (> 250 K) produce low excitation tempera¬ 
tures of the bearing molecules and, therefore, low line inten¬ 
sities. This demonstrate the importance of the isotopologue lines 
to constrain the excitation conditions of molecular lines. 

The densities and temperatures derived from our two-phase 
model for CO also fit the excitation conditions for HCN and 
HCO^, but require lower Ailing factors and molecular column 
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Fig. 13. Two component fit of the line spectral energy distribution of the (from left to right) CO, HCN, and HCO^ species, for the (from top 
to bottom) spectra at offset position (-65",+31") toward the peak of the HCN 7 = 8 —» 7 intensity, (-40",+18") toward the peak of the 
*^CO y = 16 ^ 15 intensity, (-60",-30") toward the PDR, and (-130",+30") toward the molecular gas in M17 SW. The fit of the '^C bearing 
isotopologues is shown in the insets. The cold and warm components are shown in dotted and dashed lines, respectively. The error bars are as in 

Fig.m 


densities. While two components are clearly needed to model the 
excitation towards the selected positions, relaxing the constraint 
on the isotope ratio a single component would also fit the HCN 
and HCO^ LSED towards the HCN-peak, as well as towards the 
PDR. These alternative solutions are presented in Appendix iBl 
The HCN and HCO^ LSEDs towards the HCN-peak can be fit 
using the ambient conditions of the warm component of the CO 


LSEDs, and using the same isotope ratio of 50. However, a lower 
ratio of 30 and 20 would lead to a better fit of the H^^CN and the 
H'^CO^ lines, respectively. The LSEDs towards the PDR, in¬ 
stead, can be fit using only the ambient conditions of the cold 
component of the CO, but it requires a higher isotope ratio of 75 
to fit both the H'^CN and H'^CO^ lines. 
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Table 4. LSED fit parameters” for the spectra toward the HCN-peak Table 7. LSED fit parameters” for the spectra toward the M17-W posi- 
position (-65",+31"). tion (-130",+30"). 


Parameter 

CO 

HCN 

HCO+ 

Parameter 

CO 

HCN 

HCO+ 


1.00 ±0.09 

0.40 ± 0.05 

0.40 ± 0.05 

<I>™h(‘"C) 

1.00 ±0.00 

0.80 ±0.06 

0.80 ± 0.09 

UroldQii) [cm“^] 

4.50 ± 0.47 

4.50 ± 0.47 

4.50 ± 0.46 

"coMiHz) [cm-q 

4.10 ±0.00 

4.10 ±0.35 

4.10 ±0.49 

Toold [K] 

40.00 ± 3.78 

40.00 ± 5.37 

40.00 ± 4.83 

Tcou [K] 

60.00 ± 0.00 

60.00 ±5.17 

60.00 ± 7.12 

Ncotd [cm--] 

18.90 ± 1.12 

16.60 ± 1.06 

16.40 ± 1.08 

Ncold [cm-^] 

19.00 ± 0.00 

16.20 ±0.53 

15.40 ± 1.01 


0.35 ± 0.04 

0.15 ± 0.01 

0.15 + 0.02 


0.40 ± 0.00 

0.20 ± 0.02 

0.20 ± 0.02 

n,™rai(H2) [cm-^] 

6.00 ± 0.73 

6.00 ±0.52 

6.00 ± 0.48 


5.80 ±0.00 

5.80 ±0.39 

5.80 ±0.54 

r„.„™ [K] 

130.00 ± 11.59 

130.00 ± 15.33 

130.00 ± 14.03 

r,,™,. [K] 

80.00 ± 0.00 

80.00 ± 8.61 

80.00 ± 8.79 

^warm [cm "] 

18.40 ±0.37 

15.40 ±0.54 

15.20 ±0.81 

[cm-q 

18.90 ±0.00 

15.10 ± 1.21 

15.00 ± 0.56 


1 .00±0.11 

0.40 ± 0.05 

0.40 ± 0.03 


1.00 ±0.00 

0.80 ±0.08 

0.80 ± 0.09 


0.35 ± 0.04 

0.15 ± 0.02 

0.15 ±0.01 


0.40 ± 0.00 

0.20 ± 0.02 

0.20 ± 0.02 

AV(‘2C) [kms-1] 

4.60 

7.50 

6.00 

AVi'-C) [kms-1] 

6.90 

6.10 

7.10 

AV(‘^C) [kms-1] 

3.50 

3.90 

3.70 

Ay('^C) [kms-1] 

5.50 

6.10 

4.40 


Notes. The density and column density values are given in logiQ scale. Notes. The density and column density values are given in logio scale. 


Table 5. LSED fit parameters” for the spectra toward the CO-peak po- Table 8. Mass and size associated with the cold and warm *^CO com- 
sition (-40",+18"). ponents. 


Parameter 

CO 

HCN 

HCO+ 


0.50 ± 0.05 

0.40 ± 0.05 

0.40 ± 0.04 

Hcm(H2) [cm-^] 

4.80 ±0.61 

4.80 ±0.46 

4.80 ±0.58 

Tcold [K] 

90.00 ± 12.59 

90.00 ±8.21 

90.00 ± 10.18 

Afcold [cm-2] 

18.80 ±0.35 

15.30 + 0.95 

14.40 ± 1.33 


0.10 ±0.01 

0.10 ±0.01 

0.10 ±0.01 


5.70 ±0.35 

5.70 ±0.62 

5.70 ±0.50 

[K] 

240.00 ± 29.77 

240.00 ± 26.55 

240.00 ± 26.54 

[cm-^] 

18.50 ± 1.23 

15.30± 1.10 

15.10 ±0.87 


0.50 ± 0.06 

0.40 ± 0.04 

0.40 ± 0.04 


0.10 ±0.01 

0.10 ±0.01 

0.10 ±0.01 

AV(‘-C) [kms-1] 

4.60 

4.60 

4.30 

AV(‘^C) [kms-1] 

3.90 

3.60 

3.20 


Notes. The density and column density values are given in togm scale. 

Table 6. LSED fit parameters" for the spectra toward the southern PDR 
position (-60",-30"). 


Parameter 

CO 

HCN 

HCO+ 


0.90 ± 0.07 

0.35 ± 0.04 

0.35 + 0.04 

ntoM(H2) [cm-^] 

5.30 ±0.36 

5.30 ±0.49 

5.30 + 0.62 

[K] 

60.00 ±5.32 

60.00 ± 6.82 

60.00 ± 7.09 

A/eM [cm-^] 

18.95 ±0.47 

15.40 + 0.85 

14.80 ± 1.13 


0.20 ±0.01 

0.15 ±0.02 

0.15 ±0.02 

«ra™(H2) [cm-^] 

5.80 + 0.37 

5.80 ±0.65 

5.80 ±0.56 

r„.a™ [K] 

110.00 ± 11.56 

110.00 ± 12.33 

110.00 ± 13.42 

Afnann [cm-^] 

18.60 ±0.31 

15.10± 1.10 

14.70± 1.11 


0.90 ± 0.09 

0.35 ±0.04 

0.35 ± 0.04 


0.20 ± 0.02 

0.15 ±0.02 

0.15 ±0.02 

AV(‘2C) [kms-1] 

4.30 

5.00 

4.50 

AVC'^C) [kms-1] 

4.30 

3.80 

3.90 


Notes. The density and column density values are given in log\o scale. 

The need for just a cold component in the HCN and HCO^ 
LSEDs deep in the PDR may be plausible since the HCN J - 
8 —> 7 and HCO^ 7 = 9 —» 8 are fainter at this position, indi¬ 
cating the predominance of gas that is neither warm nor dense 
enough to excite these transitions. Since the gas dominating the 
excitation of HCN and HCO^ is relatively colder in the PDR, 
a larger isotope ratio in the HCN and HCO^ lines could be 




Mass 

Diameter 

Position” 

Offset 

Cold 

Warm 

Cold 

Warm 



[Mo] 

[Mol 

[pc] 

[pc] 

200" average 

(-100", 0") 

2.6x10-^ 

102.2 

2.0 

5.1x10-2 

HCN-peak 

(-65",+31") 

l.OxlO^ 

11.2 

1.0 

1.0x10-2 

CO-peak 

(-40",+18") 

4.0x10* 

4.0 

0.4 

2.6x10-2 

PDR 

(-60", -30") 

1.0x10^ 

10.1 

0.2 

2.6x10-2 

M17-W 

(-130",+30") 

1.3x102 

40.3 

3.2 

5.1x10-2 


Notes. Selected positions are in 25" beams. 


expected (e.g. lSmith & Adamslll980HLanger et al.lll984l) . Like¬ 
wise, the warm component of the CO LSED at the M17-W po¬ 
sition also needs a larger isotope ratio in order to better fit the 
upper limit of the '^CO J - 13 —> 12 flux. 

5. Implications 

Based on our accurate estimates of the densities and column den¬ 
sities, we can further characterize the physical state of the cloud 
in terms of its parameters and energy balance. 

5.1. Mass and size traced by '^CO 

The gas mass of the cloudlets associated with each component 
of the model can be estimated from the beam averaged column 
density as, 

N- 

Mi ^ ( 2 ) 

TtCO 

where Aheam is the area (in cm^) of the emitting region subtended 
by the source size (considered equivalent to the beam size of 
25"), O; and N, are the beam area filling factors and column 
densities of the two components, and the factor 1.4 multiply¬ 
ing the molecular hydrogen mass mn, accounts for helium and 
other heavy elements. We assumed a value Xcn = 8 x for 
the r*^C01/rH?1 fraction al abundance (e.g. lLrerking et al.lll982t 
IStutzki & Gustenlll990l) . although the average value found in 
warm star-forming region s may be larger (e.g., 2.7x10“^, as 
measured in NGC 2024 bv iLacv et al.lll994l) . 
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The diameter of the cloudlets (Dd) associated with each 
component of the models can be estimated based on the den¬ 
sity and column density as Dd ~ „(h^Xco ' niass and size 
estimated for each position, and from the average spectra ob¬ 
tained from the 200" region, are shown in Table [8] The molecu¬ 
lar mass obtained from the average spectra is about 2.6x10 '* Mq 
for the cold component, while the mass associated with the 
warm component is only about 100 Mq. The mass obtained 
from the cold component is about 80% larg e r than the mass of 
1.45x10* Mq obtained bv IStutzki & Giistera dl990h using C'^O 
ob servations, but is si milar (within ~15%) to the mass estimated 
by I Snell et ^ (1 19841) from CS observations. It is also about 6 
times larger than the mass w e found from [C ii] observations 
dPerez-Beaupuits et alJl201^ . The size of a cloudlet obtained 
from the average spectra is not scaled up in the 200" region. 
It corresponds to the average size of a clump associated with the 
estimated average excitation conditions. 

The cloudlets associated with the cold components of the 
*^CO LSED at the positions of the HCN-peak, CO-peak, and 
M17-W, have larg e r size s than the median 0.11 pc and 0.15 pc 
found by iHobso^ dl992h from HCN and HCO+ 7 = 3^2 
observations, respectively. They are also larger than the extent 
(~0.24 pc) corresponding to the 25" beam size used in our maps. 
This may only correspond to the depth of the gas and may indi¬ 
cate that the cloudlet may not necessarily be symmetric when 
compared with the size projected on the sky plane. In other 
words the cloudlets extend into the plane of the sky, not surpris¬ 
ing given the edge-on geometry of this PDR. Only the cloudlet 
size associated with the cold component of the PDR at offset 
(-60", -30") is similar to those found by Hobson et al., and 
smaller than the spatial scale corresponding to the beam size. 
The sizes of the warm components are between two and three 
orders of magnitude smaller, which means they correspond to 
either a thin layer around the colder and denser cloudlets asso¬ 
ciated with the cold component, or they are actually embedded 
hot cores. 

5.2. Analysis of the energetics 

The magneti c field along the line of sig ht (Bios) toward Ml7 was 
measured by iBrogan & Trolandl (1200 11) based on their VLA H i 
Zeeman and main-line OH Zeeman observations. The resolu¬ 
tions of their H i and OH maps are ~26" and ~22", respectively, 
which allow a direct comparison with our maps. Following their 
analysis we can now study the energetics of the cloudlets asso¬ 
ciated with our selected positions in M17 SW, using the results 
obtained from our two-phase models, and assuming that the ob¬ 
served magnetic field is conserved throughout the two phases of 
the structures associated with our selected positions. 

According to iMcKee et al.l (Il993h the energy balance of a 
molecular cloud in dynamic equilibrium can be described by the 
virial equation; 

Ps + 1^1 ^ Ms+M„ + 2T (3) 

where Vs is the external pressure term, YW\ is the gravitational 
energy, Ms is the magnetic energy associated with the static 
magnetic field (Bj), M„ is the magnetic wave energy produced 
by the fluctuating magnetic field component {B„), and 27” is the 
energy contribution from internal motions (or turbulence). From 
equating the s tatic m agnetic energy to the gravitational energy 
iBrogan et al.l (1 19991) derived the following critical static mag¬ 
netic field {Bs.crit) as a diagnostic of the importance of Bs in 


a cloud: 

Bs,cri,^Sx\0-^^Np [juG] (4) 

where Np is the average proton column density of the cloud, 
estimated as 2 x A(H 2 ), where A^(H 2 ) ~ N(CO)IXco from 
our two-phase models. This is equivalent to the magnetic crit- 
ical mass (M^) used by ma ny authors in the literature (cf., 
iMouschovias & SnitzeilllBT^ . We show the Bs.crk estimated for 
the four selected positions in Table 0 If the actual static mag¬ 
netic field Bs > Bs.crit then the cloud can be completely sup¬ 
ported by Bs (is magnetically subcritical) and further evolution 
of the cloud perpendicular to the field should occur mainly due to 
ambipolar diffusion. If Bs < Bs.crit, then Bs cannot fully support 
the cloud (the cloud is magnetically supercritical) and internal 
motions must supply additional support for the cloud to be sta¬ 
ble. 

The average OH Bios detected toward the northern conden¬ 
sation of H 2 O masers, where the HCN and CO peaks are found, 
is -300 + 80 yuG, while the Bios detected in the proximities 
of the PDR at offset position ( -60",-30") is -50 0 + 80 pG 
(see B-20UC1 and B-20HI in IBrogan & Trolandl l200il their 
Sects. 3.4.6, 3.4.7, and Fig. 16). There is no quoted measurement 
of the magnetic field for the M17- W posit ion (-130",-1-30"), but 
from Fig. 16 by IBrogan & Trolandl (1200 ll) we can infer a lower 
limit of Bios = -300 + 80 pG as well, from their HI Zeeman 
observations. 

From a statistical view point the measured Bios is actually 
the average over a large number of clouds/clumps with magnetic 
fields randomly oriented with respect to the line of sight. Hence, 
it is related to the total magn etic field (]B|) as Bios ~ |B|/2 (for a 
more detailed description see ICrutcherll 1 999t Sect. 5). Thus, as¬ 
suming that the total magnetic field is the static magnetic field, 
we would have that Bs ~ -600 + 160 pG toward the north¬ 
ern condensation (including both the CO- and HCN-peaks) and 
the M17-W position, and about -1000+160 pG toward the PDR. 
These values are larger than the Bs.crit estimated for the warm 
component at all selected positions, except for the M17-W posi¬ 
tion. The cloudlets for which Bs > Bs.crit can be fully supported 
from collapse against self-gravity (if no other forces are at play). 
The estimated total magnetic fields are comparable to the criti¬ 
cal static magnetic field estimated for the cold component of the 
CO-peak and the PDR, but the Bs in the northern condensation 
is about 30% and 25% smaller than the Bs.crit estimated for the 
HCN-peak and the M17-W position, respectively. This means 
the structures associated with the PDR and the CO-peak can be 
magnetically supported, but the magnetic field cannot support 
the cloudlets at the HCN-peak and M17-W positions. In this 
case the evolution of the cloudlets at these positions should be 
controlled in part by processes that create and dissipate internal 
motions. Note that if |B| is aligned preferentially along the LOS, 
then we would have lower fields by a factor two. 

5.2.1. External ancd internal pressure 

The interactions of the forces acting on the cloudlets can be 
compared in terms of the pressure they exert. We first esti¬ 
mate the external pressure produced mainly by the radiation 
field from the most massive O and early B ionizi ng stars in 
the NG C 6618 cluster. Based on radio observations iFelli et al.l 
(Il984h derived a total luminosity by number of H-ionizing pho¬ 
tons 2(H) = 2.9 X 10^° s ”* considering a distan ce to M17 of 
2.2 kpc. On the other hand. iPellegrini et al.l (l2007l their Table 1) 
estimated a 2(H) = 1.5 x 10^° s“* for their adopted distance 
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Table 9. Critical static magnetic field and pressure terms. 



Bs,crit 

Rh® 

n c 

^ star 

Pr 

Turbulent 

pressure'' 

Thermal 

pressure 

Total internal pressure" 

Position" 

Cold Warm 




Cold 

Warm 

Cold 

Warm 

Cold 

Warm 


[pG] 

[pG] 

[pc] 

[K cm“^] 

[K cm“^] 

[K cm-3] 

[K cm-’] 

[K cm-’] 

[K cm-’] 

[K cm-’] 

[K cm-’] 

HCN-peak 

993 

314 

1.06 

1.0x10^ 

9.1x10" 

1.1x10’ 

3.3x10’ 

1.3x10" 

1.3x10’ 

2.1x10’ 

4.7x10’ 

CO-peak 

789 

395 

0.91 

1.4x10'' 

9.1x10" 

2.1x10’ 

1.6x10’ 

5.7x10" 

1.2x10’ 

3.6x10’ 

2.9x10’ 

PDR 

1114 

498 

1.34 

6.5xl0« 

9.1x10" 

5.8x10’ 

1.8x10* 

1.2x10’ 

6.9x10’ 

7.9x10’ 

2.6x10’ 

M17-W 

1250 

993 

1.58 

4.6x10" 

2.5x10" 

9.4x10" 

4.7x10* 

7.6x10’ 

5.0x10’ 

1.3x10’ 

5.2x10’ 


Notes. Selected positions are in 25" beams. Sky projected distance of the selected positions from the cluster of ionizing stars at about offset 
position (+35", +75"). To obtain the pressure terms in units of energy density erg cm“^ you must multiply them by the Boltzmann’s constant. 

The turbulent pressure is calculated from the CO line width corrected for thermal broadening of the CO molecule. The total internal pressure 
is the sum of the internal radiation pressure the turbulent pressure Pjvt, and the thermal pressure F,/,. 


of 1.6 kpc , from the luminosity of ionizing photons per star re¬ 
ported by iHanson et al.l (Il997h . For the following we have ad¬ 
justed these values to the more recent m easurement of 1.98 kpc 
for the distance to M17 (IXu et al.l 1201 Ih that we consider more 
accurate. A linear inte rpola tion betwe e n the values reported by 
iPellegrini et al.l (l2007h and iFelli et all (1 1984) give us the num¬ 
ber of ionizing photons emitted by the star cluster per second 
e(H“) = 2.39 X 10^° s-' at a distance of 1.98 kpc. 

Following the analysis bv IPellegrini et al.l (120071 their Eq.3), 
and assuming that most of the momentum in the stellar radiation 
is in the ionizing photons, we can estimate the pressure from 
the stars as Fstar - Q{H^){hv) jAnRl^ck, where {hv) is the mean 
photon energy of an O star, which we assume to be ~15 eV as 
in Pellegrini et al., Ph the distance from the star cluster to the 
position in the cloud where we want to estimate the radiation 
pressure, c is the speed of light, and k is Boltzmann’s constant. 
For the position of the star cluster we assume that of the CENla 
and CENlb stars (see Sect.[T]), at an offset of about (+35",+75") 
(cf.,Eig.l4|. The sky projected distances from this stars to the se¬ 
lected positions in Ml7 SW, as well as the corresponding stellar 
radiation pressure Pstar, are shown in Tabled 

Assuming that the (star-forming) clumps are opaque to their 
own radiation, the internal radiation pressure (Pr) due to the dust 
continuum emission (assuming the dust is well mixed in the two 
gas components) can be estimated as Pr - Tioo^mO'T^/ck, where 
TiQOfjin the continuum optical depth at 100 /J.m, cr is the Stefan- 
Boltzmann’s constant, and 7^ is the dust tempera ture. Erom the 
maps o f dust temperature and optical depth from iMeixner et aP 
(119921 their Eig. 6), we adopt = 50 K and Tioofim - 0.106 
for the HCN-peak, CO-peak and PDR positions, and = 40 K 
and Tioo/jm = 0.071 at the M17-W position. Erom the values pre¬ 
sented in Table 13 we find that the internal radiation pressure Pr 
is lower than the external stellar pressure P^tar at all the selected 
positions, except at the PDR. If these were the only forces acting 
on the cloudlets, they would collapse due to the external radia¬ 
tion pressure and self-gravity, unless the magnetic field and other 
forces support them. 

Other sources of internal pressure important to take into ac¬ 
count are the internal motions. We consider two sources of in¬ 
ternal motions: turbulence and the kinetic temperature. These 
two sources of motions are responsible for the CO line widths 
in individual cloudlets. However, it is well known from the lit¬ 
erature that turbulence (that can have several origins) dominates 
the broadening of the CO lines in warm molecular clouds. In 
our case, the average line width (EWHM) of the cloudlets we 
observe at the four selected positions is 4.5 km s ', of which 
the highest kinetic temperature of 240 K, obtained for the warm 
component of the CO LSED observed at the position of the 


HCN-peak, contributes with only ~0.63 kms^' (i.e., about 14%) 
when considering the thermal velocity contributing to the line 
widths. 

The pressure exerted by the thermal and turbulent mo¬ 
tions can be estimated as P,/, - n(H 2 )T K cm“^ and P/yr = 
n(H 2 )mH 2 (cr^j-)/k K cm“^, where 7 is the kinetic temperature, 
k is the Boltzmann’s constant, and ct/yt is the non-thermal com¬ 
ponent of the one dimensional velocity dispersion of an indi¬ 
vidual cloudlet. Other authors use the total velocity dispersion 
(assumed eq ual to the li ne width) to estimate the total internal 
pressure (cf.. lBlitjll991l) . However, we prefer to separate these 
terms, since many other authors use the thermal pressure (alone) 
in their analysis, and to compare their relative contributions in 
the two components of our models. The estimates of these terms 
are presented in Table |3 Eor the selected positions, P^t IP th 
ranges between 4 and 12 for the cold component, and between 
1.4 and 9 in the warm component. Whether thermal or turbu¬ 
lent, is clear that the internal pressure in the warm component is 
higher than in the cold component. This is also reflected in the 
total internal pressure estimated for the two phases (cf., Table|3). 
Note as well that the turbulent pressure is comparable (within a 
factor two) to the external stellar radiation pressure Pstar toward 
all selected positions, except at the southern PDR. 

At first look, the total internal pressure of these cloudlets 
seem quite high. In fact, it is about four orders of magnitude 
higher than the median pressure ~3000 Kcm“^ estimated to¬ 
ward a number of molecular clouds with densities between 
a few hundreds and 10^ cm“^, and temperatures of ~50 K 


a tev v nunareas ana lU cm ana temperatures ot ~5U Jk 

(e.g., Crutchedll99 ^ .lenkins & TrinnlbOOlt ^Ifire et al.l l2003t 


iHeiles & Troland 12004 , l2005t Troland & Crutcher l2008h . Most 
of these studies were done based on observations of H i in ab¬ 
sorption and other tracers of the diffuse gas. However, higher gas 
densities of the order of 10^ cm“^ and 10® cm“^, but with lower 
temperatures (T* < 20 K) have also been inferred from CN Zee- 
man observations (ICrutcher et al.l 1 19961 Il999t lEalgarone et al.l 
l2008h . These higher densities and temperatures lead to an upper 
limit for the thermal pressure of Pth < 2 x 10^ K cm“^, which is 
similar to the values we found in M17 SW. It is important then 
to explore whether the high internal pressure estimated toward 
our selected positions can be supported by the ambient magnetic 
field and by self-gravity. 

The magnetic pressure can be estimated from the magnetic 
ener gy density as P m Ik - |Bp/(87rk) K cm“^, where 
(see ICrutcherll 1 999l Sect. 5, for a detailed description of this ap¬ 
proximation). Using the same average values of Bios as above, 
we find that the magnetic pressure is ~8xl0^ K cm“^ in the 


the proximities of the PDR position. Thus, the magnetic pres- 
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Table 10. Sonic and Alfenic Mach numbers. 


Position 

m. 

Cold 

r 

Warm 

Ma 

Cold Warm 

Cold 

P 

Warm 

HCN-peak 

8.991 

4.987 

0.770 

4.327 

0.015 

1.506 

CO-peak 

5.994 

3.671 

0.652 

1.838 

0.024 

0.502 

PDR 

6.862 

5.068 

1.807 

3.213 

0.139 

0.804 

M17-W 

11.012 

9.536 

0.728 

5.156 

0.009 

0.585 


sure is larger than the total internal pressure of the cold com¬ 
ponent in all the selected positions, but it would be lower than 
the internal pressure of the warm components by factors ~6, ~4, 
~1.3 and ~2.6 toward the HCN-peak, the CO-peak, the PDR, 
and the M17-W position, respectively. This mean that, from the 
pressure balance point of view, the cold components are magnet¬ 
ically supported in all the positions, but the warm components 
are not, if the magnetic held remains the same throughout all the 
gas phases. 

This seems to contradict the results we found above based on 
the critical static magnetic held criteria, where two of the warm 
components were classihed as subcritical and fully supported 
by the magnetic held. This likely means that the assumption of 
equipartition between the static magnetic energy and th e grav- 
itational energy used to estimate Bs^crit by [Brogan et al.l (Il999h 
does not hold true in our case. In fact, the equipartition assump¬ 
tion YW\ - A4s does not take into account the effects of the 
thermal and turbulent pressure, which indeed play a signihcant 
role in the energy balance of dense and warm gas, as we have just 
shown. Since the turbulent pressure dominates the total internal 
pressure at all the selected positions, it is relevant to investigate 
what causes non-thermal turbulence in the cloudlets. 

5.2.2. Origin of the turbulence 

We explore two possible mechanisms that can cause non-thermal 
motions. First, the presence of outhows. The HCN lines (and in 
a lower degree the HCO^), show asymmetries in the line prohles 
that are not present in the H^^CN and H'^CO^ isotope lines, at 
least at the S/N level of the isotope lines. In particular the J„p > 3 
lines at the HCN-peak (top panel in Figs.fTTIand fTSIl show a top- 
flat shape, while most of the HCN and HCO^ lines in the other 
positions show an asymmetric line prohle with a red-shifted peak 
stronger than the blue-shifted emission. At hrst look this line 
prohles are characteristic of outhowing motions. However, we 
know the HCN lines are optically thick and can be affected by 
self-absorption. Besides, due to the highly clumpy structure of 
Ml7 SW the asymmetric line prohles could be due to different 
cloudlets overlapping in velocity along the line of sight. How¬ 
ever, the current S/N of the H'^CN and H^^CO^ lines is not high 
enough to conhrm this. Outhowing motions, on the other hand, 
need to be conhrmed by detection of additional molecular tracers 
of this process (e.g., SiO, SO, CH 3 CN, HNCO) at the selected 
positions. 

Due to the relatively strong magnetic held measured toward 
M17 SW, a second possible source of turbulence is Alfvenic 
waves or MHD waves generated in a turbulent velocity held. To 
explore this alternative we estimate and compare the sonic and 
Alfvenic numbers. The sonic Mach number is - yficrlcs, 
where cr is the observed one-dimensional velocity dispersion 
(cr = Ay/ and AV the FWHM line width), and Cj = 

{kTtlpY^'^ is the isothermal sound speed (with p - 2.33mH the 
mean particle mass for a region of molecular hydrogen includ¬ 


ing 10% helium). The Alfvenic Mach number is niA — '^ctIVa, 
where Va — |B|/( 47 rp)'^^ is the Alfven speed, |Z?| s; 2Bios the 
strength of the magnetic held, and p = 1 . 4 mH,n(H 2 ) is the gas 
density accounting for helium and other heavy elements. From 
these Mach numbers we can also estimate the ratio of thermal 
to magnetic pressures, JSp - 2 (myi/mj)^ = 2{cslVAf', which is 
a crucial parameter in theory or supercomputer simulations of 
the structure and evolution of magnetic clouds. The estimates of 
these parameters for the cold and warm components are shown 
in Table [TOl 

The internal motions at all the selected positions are super¬ 
sonic {iris > 1) by factors between 4 and 10 in the warm com¬ 
ponents, and by factors between 6 and 11 in the cold compo¬ 
nents. Instead the magnetic waves are sub-Alfvenic {niA < 1) 
in the cold components toward all the positions, except at the 
PDR where it is super-Alfvenic by about factor two. The warm 
components are all super-Alfvenic by factors between 2 and 5. 
The Mach numbers obtained for the c old components are similar 
to those obtained by ICrutcl^ (1 1 999h toward M17 SW and other 
molecular clouds. Supersonic motions would be attributed to hy¬ 
drodynamic turbulence, but Alfvenic Mach numbers niA < 1 in 
the cold components of three positions, indicates that motions 
are Alfven waves (or MHD turbulence). In fact, the plasma pa¬ 
rameter /3p < 1 in both cold and warm components toward most 
of the positions indicates that magnetic pressure dominates ther¬ 
mal pressure. Only toward the HCN-peak seems the warm com¬ 
ponent dominated by thermal pressure. These results agree with 
the previous estimates of the magnetic pressure Pm/^ obtained 
from the average magnetic held measured toward the selected 
positions compared to the explicitly estimated thermal pressure 
Pth (see Sect. l5.2.Tl and Table|9]l. 

What remains to be determined then is whether the net con¬ 
tributions of Ms, Mw and T can support the cloudlets at the 
selected positions against gravity 'W. 

5.2.3. Energy balance in the cloudlets 

The virial terms in Eq.® are 7“ = (3/2)Mcr^, Ms - 
(l/3)h|BpP^, andTF = (3/5)flGM^/P, for the kinetic, magnetic, 
and gravitational energy, respectively. In this case cr is the total 
(i.e., thermal and non-thermal) one-dim ensional velocity disper¬ 
sion (cr = Avfwhm/ V87n(2)). Following lCrutcheil(ll999l) we use 
the values a = 1.2 and h - 0.3 to correct the virial terms for 
the geometry, central density concentration, and magnetic gradi¬ 
ent effects (see also Tomisaka et al. 1998 for a detailed estimate 
of these parameters). Because an equipartition between internal 
motions and is likely to be at play in M17 SW we also as¬ 
sume that t he magnetic wave e nergy is half the kinetic energy, 
Mv, - T (iBrogan et alJll99^ . Having all these terms of the 
virial equation allow us to estimate the term of the external pres¬ 
sure from Eq.Q, assuming dynamical equilibrium. The values 
estimated for all these terms are shown in Table [TT] 

The static magnetic energy dominates all the energy terms of 
the cold components in most of the selected positions, while the 
gravitational energy is between one and two orders of magnitude 
lower than Ms ■ The static magnetic energy of the warm compo¬ 
nent is one or two orders of magnitude lower than 'W, mean¬ 
ing that the wave magnetic energy and internal motions support 
the cloudlets associated with the warm components. Only at the 
PDR position the gravity in the cold component dominates (by 
factor ~7) the static magnetic energy. The kinetic energy of the 
cold component at the PDR is ~40% larger than the static mag¬ 
netic energy, indicative of the internal motions playing the mayor 
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Table 11. Energy terms of the Virial equation. 



Ps 

Ms 

M 


2T 

|1T| 

Position" 

Cold 

Warm 

Cold 

Warm 

Cold 

Warm 

Cold 

Warm 

Cold 

Warm 


[ergs] 

[ergs] 

[ergs] 

[ergs] 

[ergs] 

[ergs] 

[ergs] 

[ergs] 

[ergs] 

[ergs] 

HCN-peak 

1.3x10'*'^ 

2.3x10^^ 

1.0x10^^ 

1.0x10'^* 

1.1x10'^® 

1.3x10*® 

2.3x10*® 

2.5x10*® 

1.2x10*® 

1.5x10*® 

CO-peak 

2.0x10^'’ 

1.3xl0'*^ 

6.6xl0'*^ 

1.7x10'*2 

4.6xl0'*^ 

4.6x10** 

9.1x10*® 

9.1x10** 

4.9x10** 

7.8x10*® 

PDR 

2.4x10^'’ 

2.5x10'*^ 

l.bxlO'*^ 

4.6x10'*^ 

1.0x10'*® 

1.0x10*® 

2.0x10*® 

2.0x10*® 

7.1x10*® 

4.9x10** 

M17-W 

3.4x10^* 


3.3X10'** 

1.3x10'*3 

3.2x10'*® 

1.0x10*® 

6.5x10*® 

2.1x10*® 

6.2x10** 

3.9x10*® 


Notes. Selected positions are in 25" beams. 


role supporting the gas against gravity at this position, despite 
the Bios magnetic field being the strongest close to the PDR. 

Since we are using an average magnetic field measured to¬ 
wards the northern condensation of H 2 O masers and close to the 
PDR, Bs may actually be stronger (or weaker) toward some of 
the selected positions. We explore scenarios where the magnetic 
held can be weaker or stronger than considered here. 

5.2.4. Alternative scenarios 

One possibility that the actual magnetic held is weaker than we 
have estimated, is if the held lines are mostly aligned along the 
line of sight and just a relatively weak held is parallel to the 
plane of the sky. If this holds true the total magnetic held would 
(approximately) be the observed Bios, and the magnetic pressure 
and the static magnetic energy would be a factor 3 lower than 
estimated. In this case the magnetic energy would still be suffi¬ 
cient to support the gas in the cold components against collapse 
for most of the selected positions, except for that in the southern 
PDR. 

On the other hand, there are two arguments favoring stronger 
magnetic helds than we estimated. First, there is evidence for 
increasing strength (between factor 2 and 4) in the magnetic 
held in small scale structures, as observed with higher reso¬ 
lution (22") with respect to ear lier Zeeman observation s done 
with coarser resolutions (60") (iBroean & Trolandll200lh . Sec¬ 
ond, and related with the previous argument, the maximum mag¬ 
netic held strength in a cloud is expected to scale with den¬ 
sity as B ^ Boinlno)^^^ for densities n > 300 cm~^, where 
no = 300 cm“^ fo r diffuse clouds and B g is the magnetic held 
at those densities (ICrutcher et alJ 1201^ . In our case we have 
that the OH observations toward the northern H 2 O condensation 
arise from gas with densities ~ 4 x 10"^ cm“^ (iBrogan & Trolandl 
l2001h . Hence, assuming Bq - Bios - -300 fiG and no the density 
for the OH bearing gas, the magnetic held strength in the cold 
component of the, e.g., CO-peak, with density n ~ 6.3x10"^ cm“^ 
(cf.. Table ID, would he B ^ -400 fiG. The total magnetic held 
at this position would then have a strength of about -800 fiG 
(assuming B 2Bios as above). Such a strong magnetic held 
would support the cloud structures in the cold component of the 
CO-peak against self-gravity (comparing with the critical static 
magnetic held Bs^crit), but not at the HCN-peak nor at the M17- 
W positions. The magnetic held would also support the cloudlets 
of the cold component at the PDR position, as well as the warm 
components at all selected positions, since their densities are 
larger. The pressure exerted by the scaled up magnetic helds 
would match even the high total internal pressure of the warm 
component (cf., Table|9l) at all selected positions. 

Higher resolution observations, as well as high sensitivities 
to detect low Stokes I and V intensities in molecular line Zee¬ 


man observations are needed to test the hypothetically strong 
magnetic helds in dense and small structures within Ml7 SW. 

6. Conclusions 

We presented new velocity-resolved maps (2.0 pc x 3.1 pc) of 
the J - 11 —> 10, 7 = 12—» 11 and 7 = 16 —> 15 tran¬ 
sitions of '^CO, as well as the 7 = 13 —> 12 isotopo- 
logue line, observed toward M17 SW using the dual-color Ter¬ 
ahertz receiver GREAT on board of SOFIA. We also completed 
extended maps of the high density tracers HCN and HCO^ us¬ 
ing the EMIR receivers on the IRAM 30m telescope, and the 
FLASH+ and CHAMP+ receivers on APEX. 

Combining this new data set we obtained the most complete 
large-scale ^^CO, HCN and HCO^ line spectral energy distri¬ 
butions (LSEDs) reported so far on M17 SW and on any other 
star-forming region, in terms of transitions covered. 

6.1. Line intensity ratios 

We found that the HCO^/HCN 7 = 1 —» 0 line ratio is lower 
than unity in all the region mapped, while the ratios between the 
higher-7 lines is larger than unity. The ratios lower than unity 
in the 7 = 1 —> 0 cold be explained with infrared pumping of 
the low-7 transition of HCN, but IR observations are needed to 
check for absorption at 14.0 fim wavelength of the first bend¬ 
ing state of HCN. The ratios larger than unity in the higher-7 
lines cannot be explained with a standard PDR model. We con¬ 
clude that the most likely mechanism that can produce the bright 
HCO^ emission is a lower excitation temperature of HCN, which 
requires larger critical densities than HCO^. 

6.2. Excitation conditions from LSEDs 

The excitation conditions of four selected positions in the maps 
were determined by fitting a two-phase non-LTE radiative trans¬ 
fer model to the CO, HCN and HCO^ LSEDs. The different 
LSED shapes associated with the selected positions are indica¬ 
tive of the distinct excitation conditions that dominate the warm 
component of the models needed to fit the high-7 CO lines ob¬ 
served with SOFIA/GREAT. The gas temperature exciting the 
gas in the cloudlet associated with the peak emission of the 
'^CO 7 = 16 ^ 15 line is ~240 K, in agreement with previ¬ 
ous estimates from the literature. This high temperature is also in 
line with the proximity of the M17 ultra- compact H ii region, an 
H 2 O maser, and at least three deeply embedded X-ray sources. 

The HCN and HCO^ LSEDs can be ht using the same am¬ 
bient conditions found for the CO lines, but with lower beam 
filling factors and column densities. All LSEDs could be fit us¬ 
ing the same isotope ratio of 50 for the three species, in line with 
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what would be expected for equilibrium chemistry of evolved 
gas. However, lower isotope ratios of 20 and 30 would lead to 
a better fit of the H'^CN and H'^CO^ lines. Inhomogeneities 
and non-equilibrium chemistry (in the proximity of YSOs, hot 
cores and UC-H ii regions) may explain an enhanced abundance 
of the rare isotopologues. An alternative solution using only a 
cold component, and a larger isotope ratio of 75, is also possible 
for the HCN and HCO^ lines at the southern PDR position. 

The diameter obtained for the warm components indicates 
they may correspond to a thin layer around the cold cloudlets, or 
to an embedded hot core or UC-H ii region. 

We also constrained the excitation conditions of the average 
LSEDs obtained over an area of 200 arcsec^. The density and 
temperature found for the cold and warm components of the av¬ 
erage spectra are similar to those found for the HCN-peak, but 
with a larger CO column density for the cold component. The 
shape of the CO LSEDs, however, resembles that observed at the 
offset position of the southern PDR. Therefore, the USED shape 
and ambient conditions derived for an unresolved source (as in 
the case of extragalactic sources with ambient conditions similar 
to Ml7 SW) would be dominated by a dense PDR like compo¬ 
nent, while the distinctive features of the high-7 CO lines ob¬ 
tained with SOEIA/GREAT for dense cores like the HCN-peak 
and dense and warm gas like that found in the proximities of the 
UC-H II region, would be smeared out. 

6.3. Energy balance and stability of the cloudlets 

We estimated the stability and energy balance of the cloudlets as¬ 
sociated with both components in our models, considering mea¬ 
surements of the magnetic field based on OH and HI Zeeman 
observations toward M17 SW, as well as the external radiation 
pressure exerted by the cluster of ionizing stars. 

We found that the cloudlets associated with the cold compo¬ 
nent of the models are magnetically subcritical {Ms > \'W\) and 
supervirial {ITIYW\ > 1) at all the selected positions, except 
for the southern PDR. The warm cloudlets are all supercritical 
{Ms < YW\) and also supervirial. 

Considering only gravity, the critical magnetic field, and the 
static magnetic field, we found that the cold cloudlets associ¬ 
ated with the southern PDR and CO-peak positions can be mag¬ 
netically supported, but the magnetic field cannot support the 
cloudlets at the HCN-peak and M17-W positions. Thus, the evo¬ 
lution of these cloudlets is expected to be partially controlled by 
processes that create and dissipate internal motions. 

Turbulent pressure dominates over thermal pressure in both 
cold and warm components at all selected positions. The total 
internal pressure of the cold components is between a factor of 
two and three larger than the external radiation pressure at four 
selected positions, and is one order of magnitude larger at the 
southern PDR position. 

The gas velocities of both components are supersonic at all 
selected positions by factors between 4 and 11. The velocities of 
the cold component are sub-Alfvenic at all the positions, except 
at the southern PDR. The velocities of the warm components 
are super-Alfvenic by factors between 2 and 5 at all selected 
positions. This suggests that internal motions of at least the cold 
components are due to MHD waves. The ratio of thermal to mag¬ 
netic pressures < 1 in both gas components indicates that 
magnetic effects dominate thermal effects at most of the selected 
positions. 

Assuming that the static magnetic field remains constant 
throughout all the gas phases, the magnetic pressure can sup¬ 
port the total internal pressure of the cold components in all the 


positions, but it cannot support the internal pressure of the warm 
components. 

If the magnetic field scales with density as B oc then 
the actual magnetic field associated with the cold components 
of our two phase model would be strong enough to support the 
cloudlets against gravity at two of the selected positions. The 
same holds true for the warm components at all selected posi¬ 
tions. This mean that ambipolar diffusion would dominate the 
evolution of these cloudlets. 
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Appendix A: Gaussian fit of spectrai iines 

The line profiles of the observed '^CO and '^CO lines are the 
result of several components overlapping along the line of sight 
and in the velocity space, which are affected by self-absorption, 
mainly in the lower-7 transitions. The lower- and mid-7 '^CO 
lines can be fit with several Gaussian components, but the 
higher-7 (7„p >11) can be fit with just one or two components. 
We selected the Gaussian component of the lower- and mid-7 
lines based on the central velocity and line width obtained for 
the single Gaussian used to ht the '^CO 7=16^15 line. 

With the exception of the 7=1^0 transitions, the higher- 
7 H'^CN and H’^CO^ lines seem to consist of a single compo¬ 
nent that can be ht with one Gaussian. Because the main HCN 
and HCO^ line prohles show more structure than that of their 
rare isotopologues, the single component of bearing lines 
is indicative of self absorption in the main lines, at least at the 
rms level of our H'^CN and H'^CO^ spectra. Therefore, we ht 
a single Gaussian component to the main HCN and HCO^ lines 
masking the line centers and using the less self-absorbed line 
wings, in order to reproduce the missing hux. 

The Gaussian parameters of the ht is summarized in Ta¬ 
bles [AT] |A2] and |A3] for all the observed transitions at the four 
selected positions. 
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Table A.l. Gaussian fit parameters for the spectra toward the HCN-peak Table A.2. Gaussian fit parameters for the spectra toward the CO-peak 
position (-65",+31"). position (-40",+ 18"). 


Transition 

/ T,„tdV 
[Kkms-‘] 

■j^peak 
^ mb 

[K] 

Vo 

[kms^'l 

AV (FWHM) 

[km s“‘] 

Transition 

f T,„bdV 
[Kkms-‘] 

rppeak 

mb 

[K] 

Vo 

[kms“‘] 

AV (FWHM) 
[kms“‘] 



‘^CO 






'^CO 



/= 1 ^ 0 

140.4±135.2 

33.77±3.37 

20.14± 0.50 

3.91± 1.27 

y = 1 - 

^0 

128.3±4.3 

25.99±1.29 

20.06± 0.06 

4.64± 0.17 

y = 2 1 

175.7± 10.6 

42.02±2.88 

19.65± 0.03 

3.93± 0.13 

7 = 2- 

-A 1 

126.8± 29.6 

30.99±7.93 

19.45± 0.06 

3.84± 0.40 

/ = 3 ^ 2 

289.3± 1.2 

68.02±0.37 

19.95± 0.01 

4.00± 0.01 

7 = 3- 

2 

215.1±0.3 

40.40±0.07 

20.09± 0.00 

5.00± 0.00 

J = 4~>3 

275.9±5.0 

64.09±1.38 

19.57± 0.02 

4.05± 0.05 

7 = 4- 

3 

214.9± 2.9 

38.98±0.86 

19.99± 0.03 

5.18± 0.09 

J = 6^5 

256.8± 0.0 

48.46±0.00 

20.67± 0.00 

4.98± 0.00 

7 = 6- 

->5 

327.1±0.5 

65.24±0.10 

20.60± 0.00 

4.71±0.00 

/ = 7 6 

329.1± 13.9 

60.46±2.83 

20.67± 0.03 

5.11± 0.10 

7 = 7->6 

291.2± 0.0 

63.87±0.03 

20.10± 0.00 

4.28± 0.00 

y = 11 ^ 10 

158.8± 31.5 

31.94±7.13 

20.79± 0.11 

4.67± 0.47 

7= 11 

^ 10 

202.9± 6.2 

39.29±1.81 

20.34± 0.07 

4.85± 0.17 

y = 12 ^ 11 

148.6± 8.4 

25.99±2.11 

19.98± 0.15 

5.37± 0.31 

y = 12 

^ 11 

123.2± 9.6 

25.09±2.52 

19.77± 0.17 

4.61±0.29 

7= 13 ^ 12 

139.6± 7.7 

25.72±2.15 

20.68± 0.13 

5.10± 0.32 

7= 13 

^ 12 

104.5±2.6 

20.14±0.80 

20.46± 0.06 

4.88± 0.15 

y= 16^ 15 

36.4± 1.4 

8.20±0.52 

20.42± 0.08 

4.17± 0.21 

y = 16 

^ 15 

53.8±3.0 

14.50±1.06 

19.71±0.05 

3.49± 0.16 



'^co 






■^co 



y = 1 ^ 0 

64.6± 0.7 

15.51±0.29 

19.37± 0.02 

3.92± 0.06 

y = 1 - 

^0 

24.7± 0.8 

5.82±0.32 

19.64± 0.06 

3.99± 0.17 

y = 2 ^ 1 

129.3± 29.5 

31.42±7.58 

19.60± 0.05 

3.87± 0.30 

7 = 2- 

-A 1 

109.1±5.2 

29.14±1.76 

19.90± 0.00 

3.52± 0.13 

J = 3->2 

104.9± 0.1 

26.26±0.04 

19.61±0.00 

3.75± 0.00 

7 = 3- 

-> 2 

72.0± 2.0 

18.27±0.67 

20.02± 0.04 

3.70± 0.09 

7 = 6^5 

99.9± 0.4 

18.72±0.12 

18.88± 0.01 

5.01±0.02 

7 = 6- 

-^5 

134.5± 2.2 

33.30±0.61 

20.40± 0.00 

3.79± 0.02 

7= 13 ^ 12 

8.3± 1.0 

2.55±0.44 

21.22± 0.19 

3.09± 0.39 

7= 13 

^ 12 

10.6± 1.2 

2.42±0.45 

20.48± 0.25 

4.14± 0.58 



HCN 






HCN 



7=1^0 

104.4± 2.4 

23.71±0.71 

18.91±0.05 

4.14± 0.08 

7=1- 

^0 

43.8± 0.7 

10.48±0.25 

18.83± 0.04 

3.93± 0.06 

y = 3-^2 

130.4± 1.3 

15.39±0.24 

19.24± 0.03 

7.96± 0.09 

7 = 3- 

-a2 

45.7± 0.6 

7.80±0.18 

19.14± 0.04 

5.51±0.10 

7 = 4^3 

86.0± 0.3 

11.73±0.07 

19.39± 0.01 

6.90± 0.03 

7 = 4- 

^ 3 

33.0± 0.4 

6.30±0.11 

19.26± 0.02 

4.93± 0.06 

y = 8 ^ 7 

10.7± 0.1 

2.09±0.05 

19.65± 0.04 

4.80± 0.10 

7 = 8- 

-A 7 

5.2± 0.1 

1.48±0.06 

19.11± 0.05 

3.34± 0.11 



H'^CN 






H'^CN 



7=1^0 

8.8± 0.5 

1.98±0.20 

18.14± 0.11 

4.20± 0.34 

7=1- 

-aO 

3.0± 0.5 

0.86±0.26 

18.01± 0.32 

3.27± 0.76 

J = 3->2 

11.0± 0.6 

2.54±0.25 

19.38± 0.12 

4.10± 0.31 

J = 3->2 

5.1±0.7 

1.25±0.29 

18.96± 0.27 

3.88± 0.69 

J = 4~^3 

6.3± 0.1 

1.64±0.07 

19.43± 0.04 

3.63± 0.11 

7 = 4^3 

2.5±0.1 

0.73±0.07 

18.95±0.09 

3.31±0.23 



HCO+ 






HCO+ 



7=1^0 

115.5±0.7 

28.44±0.28 

18.65± 0.01 

3.82± 0.03 

7=1- 

^0 

41.3± 1.1 

10.72±0.44 

19.00± 0.05 

3.62± 0.11 

J = 3->2 

157.7± 0.6 

22.45±0.14 

20.38± 0.01 

6.60± 0.03 

7 = 3- 

-a2 

57.7± 0.7 

12.31±0.27 

20.08± 0.03 

4.40± 0.08 

7 = 4^3 

147.0± 0.1 

21.48±0.03 

20.07± 0.01 

6.43± 0.01 

7 = 4- 

^ 3 

53.0± 0.4 

11.17±0.16 

19.99± 0.02 

4.46± 0.05 

7 = 9^8 

27.3± 0.4 

5.30±0.12 

19.60± 0.03 

4.84± 0.08 

7 = 9- 

-A 8 

14.1±0.6 

3.28±0.23 

19.35± 0.06 

4.04± 0.21 



h‘^co+ 






h'^co+ 



y= 1 ^0 

15.8±0.4 

4.86±0.22 

19.06± 0.04 

3.07± 0.11 

7=1- 

-aO 

4.4± 0.2 

1.79±0.16 

18.92± 0.06 

2.35±0.17 

7 = 3^2 

26.1±0.7 

6.57±0.30 

19.87± 0.05 

3.74± 0.13 

J = 3->2 

7.2± 1.5 

2.16±0.69 

19.48± 0.35 

3.16± 0.74 


Appendix B: Alternative solution for the LSED fit 

The HCN and HCO^ LSEDs at the position of the peak HCN J - 
8^7 emission, at offset position (-65",+31"), can be fit using 
only the ambient conditions found for the warm component of 
the '^CO LSED. The parameters of the model are presented in 
Table IB. l1 and the LSED fit is shown in Fig. IB. II Here we used 
the same isotope ratio of 50, but lower ratios of 30 and 20 would 
lead to a better fit of the H'^CN and H'^CO^ lines, respectively. 

The LSED of the HCN and HCO^ lines observed deeper into 
the PDR at offset position (-60",-30") can also be fit using just 
the temperature and density of the cold component found for the 
'^CO LSED, but it requires a larger isotope ratio of 75 in order 
to fit the H^^CN and H'^CO^ lines. The parameters of the model 
are presented in Table lB.2l and the LSED fit is shown in Fig. IB. 21 
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Table A.3. Gaussian fit parameters for the spectra toward the southern Table A.4. Gaussian fit parameters for the spectra toward the M17-W 
PDR position (-60",-30")- position (-130",+30"). 


Transition 

fT,„,dV 

^peak 

Vo 

Ay (FWHM) 

Transition 

f T,„bdV 

rppeak 

Vo 

AV (FWHM) 



[Kkms-‘] 

[K] 

[kms^'l 

[kms^'l 



[Kkms-‘] 

[K] 

[kms“‘] 

[kms“‘] 




‘2CO 






'^CO 



7=1- 

^0 

131.2± 16.7 

34.84±4.67 

21.42± 0.08 

3.54+ 0.15 

7=1- 

aO 

426.2± 5.8 

44.94±1.28 

20.81± 0.09 

8.91±0.22 

j = i- 

-A 1 

195.0± 26.7 

45.90±6.55 

21.28± 0.10 

3.99± 0.16 

y = 2 - 

A 1 

476.0± 26.1 

61.97±7.58 

21.57± 0.29 

7.22± 0.79 

y = 3- 

->2 

262.3± 0.0 

68.82±0.01 

21.29± 0.00 

3.58+ 0.00 

y = 3- 

A 2 

678.3±2.1 

64.33±0.45 

21.59+ 0.03 

9.91±0.06 

J = A- 

3 

273.9± 8.1 

65.37±2.29 

20.79± 0.02 

3.94± 0.07 

y = 4- 

+ 3 

648.4± 3.6 

65.30±0.84 

21.21±0.04 

9.33± 0.11 

J = 6- 

->5 

316.3+3.6 

75.79±0.96 

21.82± 0.01 

3.92± 0.02 

y = 6- 

a5 

562.7± 1.6 

66.06±0.42 

21.87± 0.02 

8.00± 0.05 

7 = 7->6 

352.5± 9.4 

80.00±2.56 

21.45± 0.03 

4.14± 0.07 

y = 7->6 

550.4± 2.7 

64.37±0.49 

21.60± 0.02 

8.03+ 0.05 

y= 11 

^ 10 

163.9+ 2.1 

31.37±0.61 

20.83± 0.04 

4.91±0.07 

y= 11 

^ 10 

157.1±3.1 

25.05±0.80 

21.24± 0.06 

5.89± 0.15 

7= 12 

^ 11 

102.4± 3.6 

18.00±0.97 

19.68+ 0.10 

5.35±0.21 

7= 12 

^ 11 

77.9± 2.4 

13.81±0.67 

20.11±0.08 

5.30± 0.19 

y= 13 

^ 12 

76.2±9.1 

14.60±2.47 

20.46± 0.28 

4.91±0.58 

y= 13 

^ 12 

73.9± 6.0 

11.52±1.59 

20.97± 0.22 

6.03± 0.68 

7= 16 

^ 15 

11.3± 1.5 

2.51±0.48 

20.15± 0.32 

4.27± 0.58 

7= 16 

^ 15 

5.1±0.8 

1.93±0.47 

21.39± 0.23 

2.50± 0.43 




'^co 






■^co 



y= 1 - 

^0 

65.2± 1.8 

12.77±0.40 

19.60± 0.06 

4.80± 0.07 

7=1- 

aO 

112.4± 0.6 

17.20±0.16 

20.47± 0.02 

6.14± 0.04 

y = 2 - 

-> 1 

135.8+3.0 

40.13±1.17 

21 .20± 0.00 

3.18± 0.06 

y = 2 - 

A 1 

342.0± 3.2 

46.82±0.68 

20.61±0.03 

6 .86± 0.07 

y = 3-*2 

151.9± 1.2 

45.56±0.47 

21 .20± 0.00 

3.13+ 0.02 

y = 3- 

A 2 

290.9+ 9.3 

39.04±1.49 

20.76± 0.09 

7.00± 0.14 

y = 6- 

-a5 

169.8± 0.2 

45.32±0.11 

21.24± 0.00 

3.52± 0.01 

y = 6- 

+ 5 

147.8± 0.1 

33.95±0.02 

20.29± 0.00 

4.09± 0.00 

y= 13 

^ 12 

3.7+ 1.3 

1.75±0.64 

18.95+ 0.37 

2 .00± 0.00 

y= 13 

^ 12 

1.1+0.3 

0.28±0.06 

21.00± 6.30 

4.0± 1.20 




HCN 






HCN 



7=1- 

^0 

70.5± 0.8 

14.84±0.23 

19.33±0.03 

4.47± 0.05 

y= 1 - 

aO 

45.9+ 1.9 

12.30±0.62 

19.22± 0.04 

3.51±0.10 

y = 3- 

-a2 

74.0± 0.9 

11.63±0.19 

19.18± 0.03 

5.98± 0.06 

y = 3- 

a2 

111.3± 6.8 

14.49±0.97 

20.76± 0.06 

7.22± 0.20 

y = 4- 

^ 3 

43.9+ 0.3 

8.04+0.09 

19.81±0.01 

5.13+ 0.04 

y = 4- 

+ 3 

49.4± 0.7 

7.36±0.14 

21.15± 0.02 

6.31+0.07 

J 

3.0+ 0.3 

0.69±0.11 

19.23± 0.22 

4.16± 0.47 

y = 8- 

A 7 

1.5+ 0.3 

0.31±0.10 

20.25± 0.58 

4.82± 1.10 




H'^CN 






H'^CN 



7=1- 

-aO 

7.3+ 0.6 

0.99+0.13 

19.17± 0.27 

6.94± 0.72 

y= 1 - 

aO 

12 .6± 1.0 

0.95+0.12 

19.48± 0.50 

12.52± 1.21 

J = 3->2 

4.2± 0.7 

1.13±0.25 

19.46± 0.32 

3.56+ 0.55 

J = 3->2 

8.2± 1.5 

0.80±0.23 

20.85± 0.85 

9.70± 2.13 

y = 4^3 

1 .6± 0.1 

0.38±0.06 

19.87± 0.20 

3.95± 0.44 

J = 4- 

+ 3 

1.9± 0.2 

0.25±0.04 

21.33± 0.37 

7.14± 0.79 




HCO+ 






HCO+ 



7=1- 

^0 

70.2± 0.5 

15.11±0.17 

19.76± 0.02 

4.37± 0.03 

y= 1 - 

aO 

91.9+ 2.3 

11.65±0.36 

20.50± 0.04 

7.41±0.13 

y = 3- 

-a2 

79.2± 0.5 

15.96+0.13 

20.44± 0.01 

4.66± 0.02 

y = 3- 

a2 

154.7± 9.1 

25.19±1.61 

21.54± 0.04 

5.77± 0.14 

y = 4- 

^ 3 

62.7± 0.2 

12.38±0.06 

20.50± 0.01 

4.76± 0.01 

y = 4- 

+ 3 

S9.4± 0.8 

14.43±0.17 

21.60± 0.02 

5.82± 0.04 

y = 9 - 

-A 8 

6.3+ 1.2 

1.47±0.46 

19.82± 0.35 

4.05+ 0.97 

y = 9- 

A 8 

9.5± 1.8 

0.91±0.24 

21.11±0.75 

9.90± 1.83 




h‘^co+ 






h'^co+ 



y= 1 - 

-aO 

6.3+ 0.3 

1.51±0.12 

19.69± 0.11 

3.96+ 0.23 

y= 1 - 

aO 

7.2± 0.2 

1.17±0.07 

21 .11±0.11 

5.83±0.25 

y = 3 2 

8.4± 0.9 

2.04±0.31 

19.87± 0.22 

3.91±0.42 

y = 3^2 

7.0± 1.0 

1.51+0.33 

22.15±0.35 

4.39+ 0.69 


Table B.l. LSED fit parameters" for the spectra toward the HCN-peak 
position (-65",+31"). 


Parameter 

CO 

HCN 

HCO+ 

KcoUi'^C) 

0.90 ± 0.09 



rioolAiii) 

4.50 ± 0.43 



Tcold [K] 

40.00 ± 4.52 



N„id [cm“^] 

19.40 ± 1.10 




0.35 ± 0.04 

0.35 ± 0.04 

0.35 ±0.02 

[cm-3] 

6.00 ± 0.63 

6.00 ± 0.42 

6.00 ± 0.32 

Twarm [K] 

130.00 ±9.62 

130.00 ± 14.11 

130.00 ± 15.78 


18.40 ±0.33 

15.10 ±0.39 

14.90 ±0.34 


0.90 + 0.10 




0.35 ± 0.03 

0.35 ± 0.04 

0.35 ±0.03 

AV('^C) [kms-1] 

4.60 

7.50 

6.00 

AVl'^C) [kms-1] 

3.50 

3.90 

3.70 


Notes. The density and column density values are given in logio scale. 
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Fig. B.l. Two component fit of the line spectral energy distribution of the CO, HCN, and HCO^ species, for the spectra at position (-65",+31") 
toward the HCN-peak. The fit of the *^C bearing isotopologues is shown in the insets. The cold and warm components are shown in dashed and 
dotted lines, respectively. The error bars are as in Fig.|9l 



Fig. B.2. Two component fit of the line spectral energy distribution of the CO, HCN, and HCO^ species, for the spectra at position (-60",-30") 
toward the southern PDR. The fit of the *^C bearing isotopologues is shown in the insets. The cold and warm components are shown in dashed and 
dotted lines, respectively. The error bars are as in Fig.|3 


Table B.2. LSFD fit parameters” for the spectra toward the southern 
PDR position (-60",-30"). 


Parameter 

CO 

HCN 

HCO+ 


0.90 ± 0.07 

0.35 ± 0.04 

0.35 ± 0.04 


5.50 ±0.60 

5.30 ±0.51 

5.30 ±0.55 

TcoU IK] 

60.00 ±8.21 

60.00 ± 7.01 

60.00 ±6.16 

NcoU [cm“-] 

18.80 ±0.67 

16.00 ±0.56 

15.50 ±0.53 


0.20 ± 0.03 



n,,onn(H2) [cm“^] 

5.80 ±0.52 



Twarm [K] 

110.00 ± 10.18 



A',™™ 

18.60 ± 1.37 




0.90 ±0.12 

0.35 ± 0.04 

0.35 ± 0.04 


0.20 ± 0.02 



[kms-1] 

4.30 

5.00 

4.50 

AVt'^C) [kms-1] 

4.30 

3.80 

3.90 


Notes. The density and column density values are given in logio scale. 
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